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Abstract 

I review the processes that shape the evolution of protoplanetary discs around young, solar-mass stars. I first discuss 
observations of protoplanetary discs, and note in particular the constraints these observations place on models of 
disc evolution. The processes that affect the evolution of gas discs are then discussed, with the focus in particular 
on viscous accretion and photoevaporation, and recent models which combine the two. I then discuss the dynamics 
and growth of dust grains in discs, considering models of grain growth, the gas-grain interaction and planetesimal 
formation, and review recent research in this area. Lastly, I consider the so-called "transitional" discs, which are 
thought to be observed during disc dispersal. Recent observations and models of these systems are reviewed, and 
prospects for using statistical surveys to distinguish between the various proposed models are discussed. 
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In this chapter I will discuss the processes that 
shape the evolution of discs around young stars. 
Such discs form as a consequence of angular mo- 
mentum conservation during th e star formation 
process (see, e.g., t he re views by IShu et all . Il987 : 
McKee fc Ostrikeil l2007h . and much of the stellar 
mass is thought to be accreted through the sur- 
rounding disc. Moreover, such discs are presumed 
to be the sites of planet formation, and as such the 
study of disc evolution has important consequences 
for theories of both star and planet formation. 
For reasons of length, in this chapter I will limit 
myself to the discussion of discs around stars of 
approximately solar mass (the T Tauri stars), but 
I note in passing that similar physics applies to 
discs observed across the stellar mass spectrum. In 
the context of this volume, this chapter falls be- 
tween the more general discussion of disc physics 
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by Lodato, and the discussion of planet formation 
and the planet-disc interaction by Klahr. More gen- 
erally, many excellent reviews of this topic can be 
found in the Protostars & Planets series (especially 
in the two most recent volumes) , and the interested 
re ader may also wis h to consult the recent review 
bv lArmitag"3(|2007al ). 

In Section 1 I discuss observations of protoplan- 
etary discs, and the constraints these observations 
place on models of disc evolution. In Section 2 I con- 
sider the processes which affect the evolution of the 
gaseous component of protoplanetary discs, and fo- 
cus in particular on models of disc photoevapora- 
tion and viscous evolution. I then go on to consider 
the evolution of the solid component of the disc (in 
Section 3), discussing the processes of grain growth 
and planetesimal formation, and the dynamical pro- 
cesses that affect dust grains. Lastly, in Section 4 I 
discuss so-called "transitional" discs, which are be- 
lieved to be crucial to our understanding of disc 
clearing. 
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1. Observational Review 

1.1. Observations of protoplanetary discs 

Discs around young stars were first directly ob- 
served in the mid-1980s, initially in the form of 
dusty debris discs suc h as those s een around Vega 
( Aumann erallll984[ ) and /3 Pic (|Smith fc Terrild . 
19841 ). and later in the form of the gas -rich discs that 



we no w term "protoplanetary" ( Sargent fc Beckwithl . 



19871 ). Advances in telescope technology in the in- 



tervening two decades has resulted in the detection, 
both directly and indirectly, of many more discs 
around young stars, and these objects are now un- 
derstood to be commonplace. Here I summarize 
the various methods of observing protoplanetary 
discs, and the constraints these observations place 
on models of disc evolution. 

Young low-mass stars are traditionally classified 
in two different ways: by the shape of their in- 
frared spectral energy distribution (SED), usually 
measured from multi-band photometry, or by the 
strength of emission lines seen by (usually opti- 
cal) spectroscopy. The SEP cl assification scheme 
was proposed by Lada |l987h and updated by 
lAndre fc Montmerld (j 19941 ). and classifies objects 
based on the slope of their SEDs at wavelengths 
longward of 2fim. Light from the central (proto- 
)star is absorbed by dust in the circumstellar en- 
vironment and subsequently re-emitted at longer 
wavelengths, so a "redder" SED, resulting from 
more material in the circumstellar environment, is 
thought to be indicative of an earlier evolutionary 
phase. The scheme is now supported by much the- 
oretical work, and while it is recognised that the 
orientation of the source with respect to the ob- 
server can have an ef fect on the observed SED (e.g., 
Whitnev et all . l2003[) . the existence of a broad cor- 



relation between SED class and evolutionary phase 
is widely accepted. Thus Class II objects are inferred 
to be stars with surrounding discs, while Class III 
sources (with near-stellar SEDs) are assumed to be 
more evolved objects which have shed their discs. 

The second classification scheme relies on mea- 
suring the strength of emission lines such as Ha. 
Young, solar-type stars with very strong emission 
lines were observed as long ago as the 1940s, and as 
the number of known such objects increased they 
were classified as T Tauri star s (hencefor th TTs), af- 



ter their " prototype" T Tau ( Jovl 19451 . see also the 
review bv lBertoutl[l989f ). It was soon realised that 



these were young, newly-formed stars, and that the 
bright emission lines (and also the bright UV emis- 
sion characteristic of TTs) were produced by the 
accretion of material on to the stellar surface. Later, 
TTs were also discovered to be bright X-ray source s 



re.g.. lKu fc ChananL[l979l:lFeigelson fc KrisaLll98l[ ). 

but the observations which detected TTs in X- 
rays also detected a further population of sources. 
These were bright in X-rays but did not show the 
strong emission lines associated with accretion. 
However, they were spatially co-located with the 
TTs and showed similar stellar properties, and so 
were christened "w eak-lined" T Tauri stars (e.g., 
Walter et al. . 198811^ . We now understand that the 
original class of objects, now referred to as "clas- 
sical" T Tauri stars (henceforth CTTs), represent 
young stars which are accreting from a circumstel- 
lar disc, while the weak-lined T Tauri stars (WTTs) 
are similar objects which have shed their discs. 
Thus Class II objects and CTTs are inferred to 
possess discs, while Class HI objects and WTTs are 
disc-less. In modern astronomical literature these 
terms are used almost interchangeably, but while 
there is a large overlap between the two classifica- 
tion schemes there is not a one-to-one correspon- 
dence: a few CTTs show Class HI SEDs, while a 



few WTTs have Class 11 SEDs fe.g.. IStrom et al 



1993t iKaas et all . [20o3) . (For further discussion of 
disc SEDs see the chapter in this volume by Wood.) 

As mentioned above, some of the most commonly 
observed characteristics of protoplanetary discs 
are signatures of gas accretion. TTs a re known to 
have strong magnetic f i elds ( '^ IkG iBasri et all 



19921 : iJohns-Krull et al.l . Il999f) . and consequently 



it is believed that the inner edges of TT discs are 
truncated by the stellar magnetic field. Inside this 
magnetospheric radius (typically ^ 10 stellar radii) 
material is channelled along the magnetic field 
lines, and essentially free-falls on to the stel l ar sur - 
face (|Ghosh fc Lamblll978HHartmann et al.lll994l) . 
Typical infall velocities are 100-300km s~^, and 
consequently we see a shock (known as the accre- 
tion shock) where the infalling material impacts the 
stellar surface. The accretion shock typic ally has 
a temperature of 2± 10 , OOOK (jCalvet fc GuUbrina . 
1998t iGullbring ei~al\ . |2000| ). much hotter than 



the stellar photosphere, and results in an ultra- 
violet (UV) excess that can be measured either 



^ Note that I keep to the modern convention of "weak-hned" 
T Tauri stars, rather than "naked" T Tauri stars as used by 



IWalter et all l ll988l) . 
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via the veilinj^ of p hotospheric absorption lines 



([Hartigan et al 



IRAC Colours of Taurus— Auriga TTs 



or by direct measur ements 
of the UV luminosity (|Gullbring et all 1 19981 . HoOfl) . 
Observations of the UV excess essentially measure 
the accretion luminosity, and if the stellar parame- 
ters (mass, radius & temperature) are known such 
observations allow the measurement of the accre- 
tion rate on to the stellar surface: typical accretion 
rates for CTT s rang e from 10-M0 - |°Mf:j yr^^ 
(|Hartigan et all Il995l : ICullbring et~al] . \l9m . In 



addition the infalling gas produces very broad 
emission lines, which are easily distinguished from 
the much narrower emission lines produced by the 
stellar chromosphere. Observations of such lines 
also allow the measure ment of accretion rates (e.g. 



MuzeroUe et al. . 200Cl[ l and, while accretion rates 



measured in this manner are somewhat more prone 
to systematic uncertainties those derived from UV 
continuum emission (as the relationships between 
line fluxes and accretion rates are empirically cali- 
brated), they allow us to detect accretion at levels 
lower than can currently be observed in the UV 



e.g., iMuzerolle et al. ■ l2003t iMohantv et al.l . [ioosh 



Observations of accretion are important as they are 
the easiest way of detecting gas in protoplanetary 
discs, but carry the important caveat that they only 
probe the innermost regions (< 0.1 AU) of the disc. 

Perhaps the most common means of observing 
protoplanetary discs is the detection of infrared 
emission from warm dust (typically a few hundred 
K) in the inner part of the disc (typically within c± 
1 AU of the central st a r.lKenvon fc Hartmannl . 119871 : 
IChiang fc GoldreichL Il997l) . The resulting excess 
emission is most commonly observed at wavelengths 
from 2-10/^m, and is readily detecte d by phq t omet- 
ric observation s (e . g.. iKenvon fc Hartmannl . 1995 



Haisch et all l2001t iHartmann et all l2005l ). Conse- 
quently, photometric surveys of nearby star-forming 
regions are capable of detecting hundreds of proto- 
planetary discs, and large surveys of this type allow 
a rapid census of the disc population. A typical 
example of such a survey (of the nearby Taurus- 
Auriga star-forming region; iHartmann et all l2005^ 
is shown in Fig. 1: the distinction between the disc- 
bearing Class II objects and the disc-less Class III 
sources is clear. 

Gas in discs is much harder to observe. Warm gas 
in the inner disc has been observed through the de- 



^ The term "veiling" describes the the "fiUing in" of ab- 
sorption Hnes in the stellar spectrum by excess continuum 
emission, in this case from the accretion shock. 
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Fig. 1. Spitzer IRAQ 2-colour plot for TTs in Taurus- Au- 
riga (data from IHartmann et al.L [20051) . The units are mag- 
nitudes: stellar photospheres have colours of zero, while red 
excesses have positive colours. CTTs (classified by H« equiv- 
alent widths) are shown as solid circles, WTTs as open cir- 
cles, and the transition object CoKu Tau/4 (see Section 4) 
as a star. Note the clear gap between the loci of the Class 
II CTTs and Class III WTTs. 

tection of molecular lines, such as CO fundamen- 



tal lines in th e infrared ([Najita et all l2003l . see also 
the review bv lNaiita et al .1 120071) . or electronic tran- 
sitions of H2 observed in the UV ( Herczeg et all 
20041 ). However, to date such studies have generally 
been restricted to relatively small numbers of ob- 
jects, and as such these data are of limited use when 
considering disc evolution. At larger radii, where 
the disc is much colder, gas is similarly hard to de- 
tect. Studies have detected rotational lines of var- 
ious mol ecules at mrn wavelengths (most notably 
CO, e.g.. lDutrev et all [l996., .2003. ), but the bright- 
est of these lines are invariably optically thick and 
therefore probe only the surface layers of the disc. 
Velocity maps of molecular lines in the brightest 
sour ces appear consiste nt with Keplerian rotation 
(e.g., Simon et all l200(/). but these too are currently 
limited to relatively small samples. 

By contrast, cold dust (~ lOK) in the outer disc 
(> 50AU) is readily obser ved in (sub-)mi ll imetr e 
continuum emission (e.g., iBeckwith et al.l . Il990l ). 
This emission is thought to be optically thin, so 
by assuming a dust opacity we are able to con- 
vert an observed continuum flux into a total dust 
mass. Standard practice then assumes a canonical 
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dust-to-gas ratio (typically 1:100) in order to infer 
a total disc mass. Disc masses measured in this 
manner ra nge from around a J u piter mass up to 
^ O.IMq (|Andrews fc Williamd . l2005h . and while 
these observations are subject to significant system- 
atic uncertainties (notably in the assumed opac- 



Taurus— Aurigo TTs 



ity and dust-to-gas ratio, iHartmann et al. , l2006h 



they provide a valuable me ans of surveying large 



ng 

numbers of TT disc s (e.g., iBeckwith et all I199C 



Andrews fc Williamd . 2005t Eisner fc Carpenter . 
20061 ). In addition, interferometric observations of 
mm continuum emission have resolved many discs 
in nearby star-forming regions, and the typical sizes 
of isolated discs are found to be ~ 100 AU (e.g., 
Andrews fc Williamsl . lioOGh . 



1.2. Observations of disc evolution 

Having discussed the various means by which 
protoplanetary discs are observed, I now discuss 
how we use such observations to constrain models 
of disc evolution. We cannot observe the evolution 
of individual objects, as the evolutionary timescales 
are much too long, so we instead draw information 
on evolution from statistical studies of disc popula- 
tions. Timescales are generally derived through the 
use of stellar ages, obtained by fitting observations 
of TTs to p re-main-sequence stellar evol utionary 
models fe.g.. Id'Antona fc Mazzitellil . Il994l ). Stellar 
ages derived in this manner are subje ct to signif 
icant systematic uncert ainties (e.g., iTout et aT 



19991 : Baraffe et al. . 2002), but the relative ages of 
different stellar populations are usually robust. 

As mentioned above, the most commonly-used 
tool in the study of circumstellar discs is near- 
infrared photometry. Large scale surveys of young 
stellar clusters provide a good base for statistical 
studies, and these studies remain one of the richest 
sources of data for disc study. A summary of the 
work of many suc h surveys is pre sented in the oft- 



cited work of Haisch et al.l (|200ll ) , which presents 



the fraction of sources with excesses (and therefore 
discs) as a function of cluster age. At stellar ages 
of less than IMyr more than 90% of objects are 
found to show infrared excesses, and the data are 
consistent with all such stars being born with cir- 
cumstellar discs. However by ages of around lOMyr 
very few significant infrared excesses are detected, 
and thus the typical lifetimes of inne r dust di s cs are 
constrained to be a few M yr (Hai sch et al.l . 12001 : 
Siciha-Aguilar et al. . 2006ah . Moreover, (sub-)mm 
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Fig. 2. IR and sub-mm observations of TTs in Taurus-Au- 
riga. The horizontal axis shows the infrared colour, in magni- 
tudes, measured between the 2MASS Kb (2.2/^m) and Spitzer 
IRAC 8.0/im photometric bands. The vertical axis shows the 
SCUBA 850/xm flux, which is thought to correlate well with 
disc mass: circles represent detections, while triangles denote 
3-cr upper limits. As in Fig. 1, filled symbols represent CTTs, 
open symbols WTTs, and the star denotes CoKu Tau/4. All 
of the accreting CTTs show strong IR excess, and the only 
CTTs not detected by SCUBA are known to be discs viewed 
edge-on. By contrast all of the WTTs show photospheric 
colours, and only two are detected in the SCUBA sample: a 
known transition object (CoKu Tau/4) and a known binary. 
This suggests that the rapid inner disc clearing inferred from 
IR observations applie s equally to the col d outer disc. [Fig- 
ure adapted from Alexander e t al.l ll2 006bl). using data fr om 
IHartmann et al.l | |2005|) and [Andrews fc Williamsl (l2005ll .] 

continuum observations of cold dust in the outer 
disc correlate very well with observation s of in- 
ner dust discs I Andrews fc Williamd . 2005 ). which 
suggests that the disc lifetimes inferred from obser- 
vations of the inner disc are indeed representative 
of the entire disc. Observations of gas are much less 
common, but measurements of accretion rates for 
individual objects show a steady decline in accre- 
tion rate with stellar age, with t he decline occur- 
ring over timescales of a few Mvr dHartmann et al 



1998t iMuzerolle et al.l . [2OO0I : ICalvet et al.l . I2OOOI ) 



However these observations allow us to constrain 
more than just the lifetimes. In surveys of low- 
mass clusters (i.e., clusters with no massive stars) 
it has long been known that very few objects are 
found between the Class II and Class HI states 
(ISkrutskie et al.', '1990*; ' Kenvon fc Hartmannl . [l995 : 
iHartmann et al., , 200 5). If we look at the infrared 
colour-colour diagrams of such clusters there is a 
striking gap between the Class II and Class HI 
loci, with few (if any) objects observed in the in- 
termediate region (see Fig. 1). Observations in the 
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mid-infrared show a similar pcap (|Persi et al 



Bontemps etall . l200lt IPadgett et all . l2006l ). and 



2OOOI : 



this suggests that the transition between the CTT 
and WTT states is extremely rapid. This conclu- 
sion is drawn from the fact that very few objects 
are "caught in the act" of clearing their discs, and 
the observations constrain the dispersal time to 
be 1- 2 orders of magnitude shorter than the life - 
time ( Simon fc Prat o[ ll99,4IWolk fc Waited . [T996h . 
This behaviour is clearly not consistent with the 
steady decline predicte d by ac cr etion disc theor; 



(iLvnden-Bell fc: Pringld . 1 1974 : Hartmann et al 



i 



19981 also the chapter by Lodato in this volume); 
instead the inner dust observations show an ex- 
tremely rapid "shut-off" in the disc emission. Mo- 
rover, observations of the outer disc show a similar 
lack of object s between the CTT and WTT states 



( Duvert et al.l . 2OOOI : [ Andrews fc Williamsl . 120051 



suggesting that the entire disc is dispersed almost 
simultaneously. This is summarized in Fig. 2, which 
shows that accretion signatures, infrared excess and 
sub-mm continuum emission all disappear contem- 
peraneously. Therefore we require that the entire 
disc, from 0.1-lOOAU, is dispersed on a timescale 
of lO'^-lO^yr after a lifetime of a few Myr. Conven- 
tional accretio n disc theory cannot e xplain such a 



rapid decline ( Armitage et al. . 1999( ). and so an- 



other explanation is required. 



1.3. Summary of disc observations 



stood. It is thought that disc evolution will result 
in a significant spread in observed properties (e.g.. 



Hartmann et al.l . Il998f ). and evolutionary models 



should be able to reproduce the range of param- 
eters discussed here. However, this cannot be the 
whole story, as we observe significant variation 
in disc parameters around stars of the same age 



in the same clusters (e.g.. iMuzerolle et al. , I2OOOI: 



Sicilia-Aguilar et al. . 2006b[ ). This suggests that the 



initial conditions for disc evolution are not uniform, 
and we note in passing that realistic models should 
also attempt to take account of the disc formation 
process (e.g., Dullemond et al. . 20061 ). 



As seen above, we are able to infer much about 
protoplanetary disc evolution from observations, 
but these observations are subject to significant 
caveats. Most important is that fact that most of 
our knowledge of disc evolution comes from obser- 
vations of dust in discs, but this dust represents 
only a small fraction of the total disc mass. This di- 
chotomy arises because the absorption and emission 
of radiation in the disc is overwhelm ingly dominated 



by th e trace dust component fe.g.. ISemenov et al 



20031 ): while gas dominates the disc mass, dust dom- 



inates the opacity. Consequently we must take care 
when comparing models to observations, as many 
disc models t reat only the ga s eous component of 



the d isc (e.g., iHartmann et al.l . Il998t IClarke et al 



20011 ). In the following sections I will discuss gas 



and dust evolution in turn, attempting always to 
anchor theoretical discussion to observations. 



Many more observations of discs have been made 
than those discussed here, but the observations dis- 
cussed here are those most relevant to disc evolu- 
tion. These observations place notable constraints 
on models of disc evolution, which are summarized 
below: 

- Disc lifetimes are typically a few Myr, with a large 
scatter. 

- CTTs and WTTs co-exist at the same age in the 
same clusters. 

- CTT disc masses range from ~ O.IM0 to < 
O.OOIMq. 

- CTT accretion rates range from ~ lO^^MQyr^^ 
to < lO-i"M0yr-i. 

- Termination of (measurable) disc accretion is con- 
temperaneous with disc clearing. 

- Disc clearing occurs rapidly (in ~ lO'^yr), across 
the entire radial extent of the disc. 

We see that all disc properties show a large scatter, 
and the origin of this scatter is not fully under- 



2. Evolution of gas discs 

2.1. Processes affecting gas disc evolution 

A number of dynamical processes affect the evo- 
lution of gas in protoplanetary discs, notably "vis- 
cous" evolutioi i due to the transport of angular mo- 
mentum (e.g., Shakura fc Sunvaeyl, 19731 : IPringle 



19811 : iBalbus fc Hawlev. . 1998, see also the chapter 



by Lodato in this volume), stellar ma gnetic pro- 



cesses such as disc w inds and jets (e.g., IShu et al 
.2000 ; Konigl fc Pudritz . 200d see also the chapter 
by Ferreira in this volum e), dynamical encounters 



with nearby stars (e .g., IClarke fc Pringld . 11993 : 
IScallv fc Clarkd . l200lh and photoevaporation of the 
disc, both by the central star and external illumina- 
tion (e.g., iHoUenba ch et all Il994 : Johnstone et al 
19981 see discussion below). The relative impor- 



tance of these processes was discussed in detail in 
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the review of HoUenbach at al. I (|200d) . which con- 
cluded that the dominant processes affecting disc 
evolution in the majority of systems are angular 
momentum transport (henceforth "viscosity") in 
the inner regions (< lOAU), and photoevaporation 
at larger radii. The other processes discussed can 
dominate over short timescales or in small regions 
of the disc (such as magnetically-launched jets very 
close to the star), or can be important for a small 
fraction of stars (such as tidal stripping by nearby 
stars, which aff ects < 10% of the stars in the Orion 
nebula cluster, IScallv fc Claris l200l[) . but do not 
dominat e the evolution of most T T discs. I refer the 
reader to lHollenbach et al 112003) for a discussion of 
the various timescales associated with these differ- 
ent processes; here I focus on photoevaporation and 
viscous evolution. I first outline the basics of photo- 
evaporation theory, and then discuss models which 
combine photoevaporation and viscous evolution. 
I conclude this Section by presenting a schematic 
model for the evolution of the gas in discs around 
(isolated) TTs, and discussing various aspects of 
current research in this area. 



2.2. Disc photoevaporation 

The term "photoevaporation" refers to mass-loss 
from a circumstellar disc due to radiative heating 
of the disc material, either from the central star 
or by an external heat source. The importance 
of photoevaporatio n was recognised as long ago 
as iBallv fc Scovilld (|l982l ). but the first models 
which treated the flow and radiative transfer prob- 



■-lOkin/s flow 



lems in a self- consistent manne r were reported by 
Hollenbach et al. (1993. ) a nd Shu et all (fl99a. an d 



then extended in detail bv lHoUenbach et al (1993). 

These models were applied with great success to 
ultracompact Hii regions around massive stars, but 
in recent years it has become clear that photoevap- 
oration is also important for low-mass stars, both in 
terms of disc evolution and planet formation. Here 
I discuss the basic physics of photoevaporation, and 
its applications to the evolution of discs around 
TTs. 

I first outline the basic principles of photoevapo- 
ration for the simplest case: that of a wind driven 
only by ionizing photons from the central star, in 
the absence of a significant stellar wind. The reasons 
for choosing this example arc twofold. Firstly, it is 
the simplest case, and almost the only one with a 
straightforward analytic solution; secondly, it is (ar- 




Fig. 3. Schematic illustration of a photoevaporative disc 
wind. Inside the gravitational radius Rg an ionized disc at- 
mosphere is formed, with scale- height H{R). Outside _Rg 
mate rial is evaporated away from the disc. (Figure adapted 
from iHoUenbach et al.lll994 12OO0I .') 



guably) the case most applicable to T Tauri stars. 
Most o f what follows rep r esent s a summary of the 
work of iHollen bach et al.l (|l994h . but for reasons of 
space and clarity I will adop t a somewhat simpler ap- 
proach. I refer the reader to lHollenbach et al.l (|l994 ) 
for a more detailed analysis. 

Consider a star of mass with a surrounding 
circumstellar disc. We assume that the disc is geo- 
metrically thin, and also that the disc is optically 
thick to Lyman continuum photons: both of these 
assumptions are valid for CTTs. Thus, if the star 
produces a strong ionizing radiation field a thin ion- 
ized layer, similar to an Hii region, is created on 
the disc surface. The ionized gas has a temperature 
of lO^K, and therefore has much greater thermal 
energy (per particle) than the cold disc material be- 
neath it. Outside some radius the thermal energy of 
the ionized gas is larger than its gravitational po- 
tential energy, and the ionized gas is unbound. As 
a result the ionized gas will flow away from the disc 
surface: this flow is known as a photoevaporative 
disc wind (see Fig. 3). 

The characteristic length scale, known as the 
gravitational radius and denoted by i?g, is found 
by equating the thermal and gravitational energies 
of the ionized gas (or alternatively by equating the 
local sound and orbital speeds), and is given by 



GM., 



8.9 



AU 



(1) 



where Cg is the sound speed of the ionized gas, typ- 
ically 10km s~^. The flow rate outside Rg depends 
on the density at the base of the ionized layer, and 
to first order the mass-loss rate per unit area is sim- 
ply evaluated as pCg at the base of the ionized re- 
gion (where p is density) . The density at the base of 
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the ionized region depends on the radial and vertical 
structure of the disc atmosphere, and consequently 
we can reduce the evaluation of the wind rate to a ra- 
diative transfer problem. In making this assumption 
we have neglected the flow structure of the wind, 
but as we will see below this makes only an order- 
of-unity difference to the calculations. Moreover, as 
we are considering only ionizing radiation the (ex- 
tremely complicated) radiative transfer problem is 
reduced to a calculation of ionization equilibrium. 

Inside Rg the scale-height of the ionized atmo- 
sphere is much larger than the scale-height of the 
(cold) disc. We can therefore approximate the struc- 
ture of the ionized layer as that of an isothermal disc 
in hydrostatic equilibrium, with scale-height H{R) 
given by 



(ciR^ 



1/2 



R<Rg (2) 



where vy^ is the Keplerian orbital velocity. The sound 
speed (temperature) of the ionized gas is constant, so 
we can substitute from Equation 1 to find a simpler 
form for the scale-height: 



H{R) = i?g 



R 

Rs 



3/2 



R < R„ 



(3) 



As mentioned above, the critical factor in deter- 
mining the flow rate is the (numb er) density at the 
base o f the ionized region, nn(i?). [llollenbach et al.l 
(jl994l ) showed that the attenuation of stellar Lyman 
continuum photons through the disc atmosphere is 
very high, and that the diffuse (recombination) field 
dominates for all radii of interesQ- If we assume 
that the diffuse flux incident on the disc surface at a 
radius R is dominated by the disc atmosphere very 
close to R (an assumption veri fied by the numerical 
work oF Hollenbach et al.ll99 4') , then we can adopt a 
Stromgren condition, equating ionizations with re- 
combinations at the base of the ionized layer: 



a^R^nliR) cx $ , R < R^ 



(4) 



In this equation a-Q — 2.6 x 10 ^■^cm'^s ^ is the Case 
B rec ombination coefficient for atomic hydrogen at 
10% (ICro3.l200rf). <& is the stellar ionizing luminos- 
ity (in units of photon s~^), and the proportionality 
accounts for the difference between the diffuse radi- 
ation field at R and the stellar ionizing flux <&. Thus 



^ Note that when hydrogen undergoes radiative recombina- 
tion approximately 1/3 of recombinations are to the ground 
state, and therefore emit another ionizing photon. 



no(i?) 



R 

Ra- 



-3/2 



R < Ra- 



(5) 



where rig is the base density at Rg. We calculate Ug 
by means of a Stromgren criterion similar to that 
above, so 



= C 



3$ 



47raBi?g 



1/2 



(6) 



where the order of unity constant C reflects the 
difference in intensity between the stellar and 
diffuse ionizing fields . The numerical work of 
Hollenbach et al.l (|l994[) fixes the value as C ~ 0.14. 



At radii beyond Rg, in the "flow" region, the ion- 
ized layer flows at approximately the sound speed. 
There is no longer a direct line-of-sight from the 
star to the top of the disc atmosphere, so the irra- 
diation beyond Rg is dominated by the diffuse field 
produced at Rg. Consequently, at radii beyond Rg 
the incident flux is simply the number of recombina- 
tions at Rg, corrected for geometric dilution. If we 
equate ionizations and recombinations at the base 
of the disc atmosphere for R > Rg, we find 



aBRlnl 

^ ^\ R 

and therefore 



R 

no{R) =ng[ — 



aei? «o(-R) > R>Rg (7) 



-5/2 



R > Rg. 



(8) 



The mass-loss per unit area from the disc, Ewind(-R), 
is given by the product of the mass density and the 
sound speed, so 

Swi„d(i?) = 2po(i?)cs = 2no{R)cs^^mn , R > Rg{9) 

where mu is the mass of one hydrogen atom, is the 
mean molecular weight of the gas (taken to be = 
1.35) and the factor of 2 accounts for mass-loss from 
both sides of the disc. There is no mass-loss within 
Rg (i.e., Swind = for i? < Rg.), and the mass-loss 
profile from the disc falls off as R~^^^ outside Rg. 
The integrated mass-loss rate is found by computing 



(10) 



Mwind = / 27ri?Ewi„d(i?)di?, 



and if we substitute the form for Swind (R) , integrate, 
and re- write in terms of parameters typical of TTs, 
we find that 
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wind 



4.4 X 10" 



1041S-1 



1/2 



(11) 



1/2 



Moyr-i. 



Thus the mass-loss profile is entirely specified by 
just 2 parameters: the stellar mass and ionizing flux. 
Moreover, the rapid fall-off in the mass-loss rate per 
unit area at radii beyond Rg means that mass-loss 
very close to Rg dominates over mass-loss at larger 
radii. 

This simplified scheme captures the qualitative 
processes well, but more recent work has resulted 
in a number of quantitative improvements to the 
model. The presence of dust in the ionized gas can 
alter the radiative transfer problem (although this 
is more applicable to mass ive stars than to TTs, 
e.g.. lRichling fc Yorkel . ll997l ). and the effects of non- 
ionizing UV radiation can also be important (e.g., 
Johnstone et al. . 19981 see also Section 2.2.1 below). 
Detailed consideration of the hydrodynamic prop- 
erties of the wind have also resulted in improve- 
ments to the wind models. The isothermal wind can 
be conside red as a Bernoulli flow ( see, e.g., the dis- 



cussion m 



Dullemond et al.l . l2007l ). and when pres- 



sure gradients are taken into we find that the "crit- 
ical radius" at which the wind is launched is in fact 



~ Rg/5 (|Liffmanl . 120031). Add i tional ly, the hydro- 
dynamic models of iFont et al. I (|2004l) showed that 
the transition between the "static" and "flow" re- 
gio ns is not as sharp as th e step-function assumed 
by ( HoUenbach et al.l . ll994l ). and find that the wind 
rate is reduced by a factor of ~ 3 (primarily due to 
the somewhat smaller, sub-sonic, launch velocity). 
As mentioned above, the qualitative behaviour of 
the wind is unchanged, but the reader s hould note 



that I will make use of the wind models of Font et al 



(|2004! ) in subsequent calculations. 



2.2.1. External irradiation: the ONC proplyds 

As mentioned above, early photoevaporation 
models tended to focu s on massive star s (e.g., 
Hollenbach et al.l . Il994l : iRichling fc Yorkel [1997), 
and only in more recent years has the application of 
these models to TTs become popular. However, it 
has long been known that photoevaporation of TT 
discs due to external irradiation can be important, 
and this is especially relevant to the case of the so- 
called "proplyds" in Orior[I]. These are a small num- 



^ The term "proplyd" was originally coined as an abbre- 
viation of PROtoPLanetarY Disc, following observations of 



ber 100) of TTs with discs observed in silhouette 
(i.e., in absorption) against the diffuse background 
emission in the Orion nebula, and are typically 
surro unded by characteristic cometary-shaped neb- 
ulae (lO'Defl et all . 119931: iMcCauehrean fc O'Delj 
19961 ). They are also bright in emission lines such 
as Ha and [Nil], and were quickly recognized as 
discs which are being photoevaporated by the UV 
emission from the nearby 0-stars a t the centre 
of the cluster. Johnstone et al. ( 19981 ) constructed 
detailed models of this process, and these models 
had great success in explaining the observed emis- 
sion. In these models the wind is primarily driven 
by non-ionizing, far- ultraviolet (FUV), photons, as 
the ionizing (EUV) photons are absorbed in the 
flow well above the disc surface. The wind rates 
are large, ~ lO~^M0yr^^ (although this depends 
strongly on the proximity of the proplyd to the 
ionizing source), and the characteristic cometary 
shape is explained by the interaction between the 
photoevaporative flow and the stellar wind from the 
nearby O-star(s). However, the combination of this 
large wind rate and the measured disc masses gave 
rise to the so-called "proplyd lifetime problem" : in 
the presence of such winds, the typical disc lifetimes 
should be so short (^ lO^yr) that we should not 
observe many, if any, proplyds. 

A neat solution to the pr oplyd lifetim e prob- 



lem has b e en pr o posed bv lAdams et al.l (|2004l ) 



and lClarkel (|2007t ). lAdams et al.1 (|2004l) noted that 



the gravitational radius for FUV-driven flows is 
typically greater than the outer disc radius (as 
FUV-heated winds have much lower temperatures, 
- 100-lOOOK, than the EUV-heated winds dis- 
cussed above). They constructed hydrodynamic 
models of the flow, and found that the wind rate 
is in fact a strong function of the disc size. If the 
disc is larger than ~ lOOAU in size then the outer 
radius is larger than the gravitational radius, and 
the wind rate is large. However, smaller discs have 
little or no material outside the gravitational radius 
and consequently produce weaker winds (launched 
radially from the outer di s c edg e) for the same ex- 
ternal irradiation. Clarke (2007) coupled the wind 
profiles of lAdams et al. (2004) to a simple viscous 
disc evolution model, and found that discs tend to 



such discs in silho uette against t he ba ckground emission in 
the Orion nebula llO'Dell et al.l. [1993'). However, more re- 
cently this terminology has become disfavoured in the liter- 
ature, and the term "proplyd" is now generally used to refer 
only to discs of this specific type: discs that are being pho- 
toevaporated by the radiation from nearby massive stars. 
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tend to evolve towards a quasi-steady state where 
the wind rate matches the viscous spreading of the 
disc. In this model small discs rapidly expand until 
the wind becomes efficient enough to limit the disc 
size, while large discs are rapidly photoevaporated 
down to smaller sizes. The dominant factor con- 
trolling disc evolution is the initial disc mass, and 



Clarkd (l2007|) find that a simple spread of initial disc 
masses can explain the observed distribution of disc 
sizes in Orion, solving the proplyd lifetime prob- 
lem. However, some uncertainties remain in this 
picture, notably that the theoretically-predi cted 



mass-loss rates (from, e.g.. lAdams et al.l . 120041) are 
still significantly smaller than those inferred from 
spectroscopic observatio ns of proplyd flows (e.g., 
Hennev fc O 'Deli Il999h . One one hand it seems 
that the low mass-loss rates predicted by models 
solve the proplyd lifetime problem; on the other 
hand, that the higher rates derived from observa- 
tions suggest that the proplyd lifetime problem still 
exists. Future studies, both observational and theo- 
retical, will hopefully resolve this discrepancy. Note 
also that, as mentioned in Section 2.1, external ir- 
radiation of this kind affects only a small fraction of 
protoplanetary discs, but that the consequences for 
those discs which are affected can be spectacular. 

2.3. Photoevaporation plus viscous evolution 

Recently, models which incorporate both vis- 
cous evolution and photoevaporatio n by the cen- 



tral s t ar have become po p ular ( e. g.. IClarke et al 



2OOII: iMatsuyama et all l2003t iRudenl . 12004 



Takeuchi et all . l2005t [Alexander et all l2006bl ) 
Models of this type show a number of attractive 
properties: here I review the basic principles of 
this c lass of models, first outlined by IClarke et al 
(l200lh . 



2.3.1. Basic concepts 

Conventional models of disc evolution ( e.g. 



Lvnden-Bell fc Pringk . 19741 : Hartmann et al.l . 



19981 see also the chapter by Lodato in this volume) 



give rise to a power-law decline in the disc surface 
density with time, and a similar decline in the mass 
accretion rate. Thus while such models are broadly 
consistent with observations of TTs, they are un- 
able to reproduce the rapid clearing observed at the 
end of the disc lifetime (see Section 1.2). The ba- 




t / yr 

Fig. 4. Decline of accretion rate (on to the star) in differ- 
ent disc models. The thi n solid line shows the e volution of 
the similarity solution of iHartmann et al.l l ll998l) . with the 
fiducial photoevaporative wind rate shown as a dotted line. 
The heavy solid line shows the accretion rate when a pho- 
toevaporative wind is incorporated into the disc evolution 
model: the behaviour follows the similarity solution at early 
times, but the wind results in a rapid termination of accre- 
tion once the disc accr etion rate falls to th e level of the wind. 
(F igure adapted from[Clarke et al.ll200ll , using the models 
of lAlexander et al.ll2006bl .') 

a disc, the accretion rate through the disc falls to 
a level comparable to the rate of mass-loss due to 
photoevaporation. The photoevaporative mass-loss 
is concentrated near the gravitational radius, so at 
this point in the evolution the outer disc, beyond 
Rg, is unable to re-supply the disc inside Rg (as 
all of the accreting material is "lost" to the wind). 
Once deprived of re-supply the inner disc de-couples 
from the outer disc, and viscosity drains the in- 
ner disc on to the star. This occurs on the viscous 
timescale of the inner disc (which is much shorter 
than the disc lifetime, see Section 2.3.3 below), and 
the inner disc therefore disappears very rapidly, in 
a manner consistent with the many of the observa- 
tions described in Section 1.2 (sec Fig. 4). To date 
this is the only class of models to have reproduced 
the observed two-timescale behaviour, and due to 
the role of the wind in precipitating this behaviour 
this model was christened the "UV-switch" . 



2.3.2. The "direct wind" 

This class of models has a number of attrac- 
tive propert i es, b ut the scenario presented by 



sic pr emise of the model proposed by IClarke et al.l 
( 200ll ) is that, at some late stage in the evolution of 



Clarke et al.l (|200lh is not valid after the inner disc 
has drained. The wind models described in Section 
2.2 assume that the disc is optically thick to ion- 
izing photons at all radii, but this assumption is 
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no longer valid once the inner disc has drained to 
a sufhciently low level. Once the inner disc drains 
the ionized atmosphere disappears, eliminating the 
diffuse radiation field, and instead we must con- 
sider the consequences of direct irradiation of the 
in ner disc edge. Thi s proces s was modelled in detail 
by [Alexander et al. ( 2006a ) ; here I summarize the 
basic results. 

Much insight into this problem can be gained from 
an analysis similar to that in Section 2.2: neglecting 
the hydrodynamic properties of the wind, and reduc- 
ing the problem to one of ionization/recombination 
balance. We again assume that the mass-loss per 
unit area from the disc at a given point is given by 
/5Cs, and again seek to evaluate the density at the 
base of the ionized layer. 

We neglect recombinations between the source 
and the disc surface, and also assume azimuthal 
symmetry (integrating over the azimuthal coordi- 
nate (f) throughout). Therefore, along any given line- 
of-sight from the source the rate of recombinations, 
A^rcc, in a volume AV^ at the ionization front must 
balance the rate of ionizing photons absorbed at the 
front, TVion- A column with polar angle 9 and an- 
gular size A9 has an area of 27rr^ sin 6 Ad, so for an 
ionizing flux $ the ionization rate at the front i J^ 



(12) 



We make the simplifying assumption that the vol- 
ume AV has a thickness equal to the disc scale- 
height H (as perturbations to the disc on scales 
shorter t han H are unlik ely to be dynamically sta- 
ble, e.g.. lLin et al.l . 19851 ) so AV is given by 



AV = 2'kRH 



rAO 
sin/3 



where (3 is the angle between the ray-path and the 
ionization front. The cylindrical radius R — rsmO, 
so the recombination rate is 

RAO 



2naBnoRH- 



sin 9 sin /3 



(14) 



We then equate the ionization and recombination 
rates and re-arrange to find 



/ $ sin^ sin /3 \ 



1/2 



(15) 



^ To avoid confusion I adopt tlie notation (r, 6, (j>) for spher- 
ical coordinates and {R, z, ip) for cylindrical; upper-case R 
represents cylindrical radius, while lower-case r denotes 
spherical radius. 



We see, therefore, that in general no{R) depends on 
the geometry of the ionization front, which in turn 
depends on the density structure of the unperturbed 
disc. This problem requires numerical solution, but 
it is instructive to frame the problem in terms of 
dimensionless scaling constants. 

We again evaluate the wind mass-loss rate per 
unit area as 



Sdiroct(-R,i) = 2nmiinQ{R,t)vi{R,t) 



(16) 



The launch velocity v\{R,t) is of order the sound 
speed of the ionized gas, and we have already made 
the simplifying assumption that at the inner disc 
edge AV has a radial thickness comparable to the 
vertical scale-height. Along the disc midplane we 
have sinO — sin/3 = 1, so the inner edge density riin 
can be expressed as 



C 



(47raB(i//i?)i„i?f„ 



1/2 



(17) 



where i?in is the radius of the inner disc edge and C is 
an order-of-unity scaling constant. We then assume 
that riQ^R) can be related to by a dimensionless 
shape function / (R/Rin), so 



no(i?) = Tlinf 



R 

R\n 



(18) 



If we write then the launch velocity as 

vi{R)^Dc,, (19) 

where D is another order-of-unity scaling constant, 
the mass-loss profile takes the form 



(13) i;dircct(^,0 = 2CD^j,miiCsnin{t)f 



R 



Rinit) 



(20) 



The total mass-loss rate is found by integrating this 
from i?in to some outer disc radius Rd, and if we 
substitute for rtjn from Equation 17 and re-scale to 
parameters typical of TTs (taking {H/ R)\a = 0.05) 
we find that 



Mdircct(i?d) = 2.35 X 10-9 Ci? 



Vio^ig 



1/2 



/ i?in \ 



1/2 



xf{x)dxMQyi-\ (21) 



We have seen in Section 1.1 that TT discs typically 
have sizes > lOOAU, so we expect that i?d ^ -Rin- 
Further, if we assume that the integral converges 
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Fig. 5. Snapshot of one of the hydrodynamic wind simula- 
tions of lAIexander et alj ll2006al ). The gas density is plotted 
as a colour scale, with the grid boundaries denoted by solid 
lines and the ionization front by a dashed line. Velocity vec- 
tors are plotted at regular intervals. 



as Rd/Rin — > oo (i.e., that the form of f{x) falls 
off faster than x"^), then we see that the mass- loss 

1 /2 

rate from the wind scales as R^^ . Consequently, we 
expect the mass-loss rate to incre ase as the inner 
edge of the disc evolves outwards. [Alexander et al 



(|2006ar ) conducted a series of hydrodynamical sim- 
ulations of this problem (see Fig. 5), and found that 
this analytic form does indeed provide a good fit 
to the wind profile. The constants CD and a were 
found to depend weakly on the disc thickess H{R): 
typical values (for H/R = 0.05) are CD = 0.235 and 
a — 2.42. Note also that the wind rate is around an 
order of magnitude larger than that of the diffuse- 
field wind, primarily due to the more efficient radia- 
tive transfer process. Consequently, we expect the 
direct wind to dominate over the diffuse wind once 
a gap has opened in the inner disc. 



2.3.3. Evolutionary timescales 

Having established the behaviour of the photoe- 
vaporative flow during the different phases of disc 
evolution, we can now consider the evolutionary 
timescales that result. [Simil ar an alyses were pre- 
sented by IClarke et al.l (|200l[ l and I Alexander et al 
(l2006bl ).] Sere we make use of the simila rity solu- 



tion presented by iHartmann et al. neglect- 
ing factors of order unity, and for simplicity adopt 
a linear viscosity law [u oc R). Consequently the 
local viscous timescale in the disc is given by 



and the characteristic scaling time for the evolution 
of the disc is 

to = ^, (23) 

where i?o is the characteristic radius of the 
similarity solution (typicall y lOAU for TTs, 
Hartmann e t al. 19.98; Clar ke et al.ll200l[ ). and v^) is 



the value of the viscosity at radius Rq. When con- 
sidering the evolution of the disc, there are three 
relevant timescales: 



(i) 



(ii) 



(iii) 



The "lifetime" of the disc (denoted by idisc): 
the time taken for the accretion rate to fall 
from its initial value to the level where photo- 
evaporation becomes significant. 
The "inner clearing" or "gap-opening" 
timescale (tinner): the time required for vis- 
cous draining of the inner (< i?g) disc after 
the photoevaporative wind creates a gap in 
disc (Clarke et al., 200ll). 
The "outer clearing" timescale (toutor): the 
time required for the (direct) wind to clear the 
outer ( > _Rg) part of the disc (Alexander et alj, 



2006bD 



The disc lifetime can be estimated simply from the 
similarity solution, and is given by 



^disc ^0 



Macc(O) 



2/3 



(24) 



wind 



where Macc(O) is the initial disc accretion rate. 
The inner clearing timescale is simply the viscous 
timescale at i?g, so 



^11" 



in 



(25) 



If we neglect viscous evolution of the disc, the time 
required to clear the outer disc to a radius R can be 
estimated as the disc mass inside R divided by the 
direct photoevaporation rate. The disc mass in turn 
is given by the product of the local viscous timescale 
and the accretion rate, so 



t„ 



Md(< R) 

Mdircct(i?) 



tu{R) 



direct 



(i?) 



(26) 



The accretion rate at which the gap opens is con- 
stant (approximately the diffuse wind rate), so the 
outer clearing timescale varies only with Mdiroct(-R)- 
We saw in Section 2.3.2 that A/dircct(^) oc i?^/^, and 
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by noting that the time required to clear the disc 
out to i?g is simply ijnncr we see that 



Evolution of surface density 



to 



tAR) ( 



^1/2 



(27) 



(i?/i?g)^^/^ < 1, so the direct wind clears the outer 
disc on a timescale shorter than the local viscous 
timescale (verifying our earlier decision to neglect 
viscosity jjj. We also see that the outer clearing time 
is longer than the inner clearing timescale by a factor 
of (i?/i?g)^/^: it takes somewhat longer to clear the 
outer disc than the inner disc. In this scenario the 
ratio of the clearing timescale to the disc lifetime is 
therefore 



touter tu{Rd) I Rd 



Rd ( Rd 



Ro \R, 



^1/2 



-1/2 



/ Macc(O) 

y M^iad 

' Macc(0) \ 



-2/3 



-2/3 



(28) 



If we select parameters typical of TTs we find that 
idisc is typically a few Myr, and the ratio of the 
clearing timescale to the disc lifetime is 1-5%. Thus 
models which combine viscous evolution with pho- 
toevaporation of the disc evolve on timescales con- 
sistent with the observational constraints discussed 
in Section 1.3. 

2.4. Schematic picture of gas disc evolution 

We now have a model for the evolution of the gas 
disc over the entire disc lifetime, and we can describe 
it as "three-stage" model for gas disc evolution (see 
Fig. 6): 



(i) 



(ii) 



Viscous phase: few Myr. In this phase, 
which encompasses the majority of the disc 
lifetime, the photoevaporative wind is neg- 
ligible and the disc evolves due to viscous 
transport of angular momentum. Most of the 
disc mass is accreted on to the star, and most 
of the angular momentum is transported to 
large radii. 

Gap-opening phase: < lO^yr. In this phase 
the (diffuse) photoevaporative wind cuts of the 
inner disc from re-supply, and the inner disc 
drains on its own, short, viscous timescale. 



^ Note that this remains true as long as MdircctC^) is an 
increasing function of R. 




Fig. 6. Snapshots of disc surface d e nsity e volution from the 
fiducial model of [Alexander et al. I l l2006bl '). colour-coded to 
identify the three stages of evolution described in the text. 
Snapshots are plotted at t = 0, 2, 4, 5.9Myr (blue), 6.0, 6.01, 
6.02, 6.03Myr (green), 6.04, 6.05. . . 6.18Myr (red). 

(iii) Clearing phase: few lO^yr. In this phase 
the inner disc is optically thin to ionizing radi- 
ation, and the disc is dispersed from the inside- 
out by the (direct) photoevaporative wind. 
This model has a number of attractive properties 
and, as discussed above, is the only model to date 
which satisfies the observed constraints on evolu- 
tionary timescales. However, many uncertainties re- 
main, and here I discuss some of the most important. 

Perhaps the biggest uncertainty in models of disc 
photoevaporation is the UV luminosity of the cen- 
tral star. The models discussed here require that the 
star produce a large flux of ionizing photons (the 
fiducial value of $ = IC^^photon s~^ corresponds to 
a luminosity of ~ 10"^°erg s~^), and moreover require 
that the ionizing luminosity is sustained through- 
out the disc lifetime. If the ionizing flux declines 
with time (as it would if it was powered by accre- 
tion) then the photoevaporative wind rate also de- 
clines, and the resulting disc evolution is inconsis- 
tent with the observational constraints discussed in 
Section 1.3 ( Matsuvama et al .. 2003; Ruden, 2004). 
However, if the ionizing luminosity is produced in 
the stellar chromosphere, and is therefore driven by 
stellar magnetic activity rather than disc accretion, 
it is expected to be approximately constant over 
the lifetime of the disc. X-ray emission from TTs 
(which is magnetically-driven) is uncorrelated with 
disc properties, and in fact some o f the most X-ray 



luminous TTs are disc-less WTTs (jFeigelson et al 



20071 ) . The solar chromosphere emits L yman contin 
uum photons at ^ lO'^^photon s~^ (e.g. jBasri et al 
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19791: lAvresl . Il997t ) and, given that TTs are much 
more magnetically active than the Sun, it is not un- 
reasonable to expect TTs to have significantly larger 
ionizing luminosities. Unfortunately direct observa- 
tions of TTs in the UV are extremely difficult to 
make, and consquently there are f ew direct estimates 
of <1> available in the literature. iKamp fc Sammar 
(|2004l ) used a scaling argument to estimate the UV- 
spectrum of a young, active G-type star, and their 
spe ctrum has an ion i zing fl ux of 2.5 x 10'''^ photon 



Alexander et al.l (|2005l ) used archival UV obser- 



vations to reconstruct chromospheric emission mea- 
sures and estimate $ for a small number of CTTs, 
and foimd <i> ~ lO^^photon s~^ to be typical. They 
also found that the Hen 1640A line, thought to be 
produced by radiative recombination, was uncorre- 
lated with any disc properties, suggesting that the 
ionizing emission is essentially unchanged during the 
TT phase. These observations, however, are limited 
to a very small sample size, and given the lack of 
observational facilities currently available in the UV 
it is unclear when better such observations will be 
possible. 

However, a great deal of indirect evidence has also 
emerged to support the h ypothesis that TT s pro- 
duce a strong ionizing flux. Font et al. (|2004^ mod- 
elled the emission-line properties of the photoevap- 
orative flow, and found that their model successful- 
Uy reproduced observed emission-line properties of 
several CTTs. More recently, observations in the in- 
frared have identified a number of emission lines that 
are thought to arise from radiative recombination in 
the disc atmosphere, most notably the [Neii] line at 
12.81/zm. This line has recently bee n detected in the 



mid-IR spectra of a number of TTs ( Pascucci et al 



20071 iLahuis et al.. . ,200 7) , and can be indicative of 



ionization of disc material by eith er EUV (HoUen- 
bach & Gorti, in prep.) or X-ray ("Glassgold et al.', 
[20071 photons. If the emission is due to EUV ioniza- 
tion, then the inferr ed ionizing fluxes are once again 



lO'^iphoton s"^ (|Pascucci et all . l2007t ): however 



if the emission is X-ray driven then it is unlikely to be 
important for disc dynamics, as X -ray heating does 
not drive a signifi cant di sc wind ('A lexander et al" " 



[2b04bf ^. Espaillat et al.l (|2007, ) found evidence for 
a correlation between [Nell] emission and stellar ac- 



® Note, however, that while the EUV fluxes are uncertain, 
the X-ray emission from these objects is readily observed. 
Consequently, any constraints on the EUV emission derived 
from [Nell] emission are technically only upper limits, as 
some fraction of the line flux is due to the (known) X-ray 
emission. 



cretion rate but no correlation between X-ray emis- 
sion and [Neil] emission, and interpreted this as ev- 
idence for (E)UV-heating of the disc atmosphere. 
However, previous theoretical models have rejected 
accretion-shock emission a s a mechanism for driv- 
ing disc photoevaporati on ( Matsuvama et al.l 120031 : 



Alexander et al.l . l2004al ). Also, the correlation covers 



only a small region of the observed parameter space, 
with a small sample size (7 detections), so its signif- 
icance remains unclear. Moreover, observations of 
disc emission lines in the mid-IR are hindered by 
the bright continuum emission from the disc, and 
consequently most of the detections of [Neil] emis- 
sion to date have been observed in so-called "transi- 
tional" discs, which are characterised by unusually 
weak continuum emission in the near- and mid-IR 
(see Section 4). It seems likely, therefore, that these 
observations are subject to significant selection bi- 
ases, and that larger samples, with better signal- 
to-noise, will be required to break the degeneracy 
between X-ray and EUV excitation. A possible so- 
lut ion to this problem h as recently been suggested 
bv iHerczeg et al. ( 2007 ). who used high-resolution 
(echelle) spectroscopy to resolve the [Neii] line in 
the nearby transitional disc TW Hya. They mea- 
sured an intrinsic line- width of ~ 20km s~^, which 
strongly suggests that the line is formed in a hot 
disc atmosphere at a radius of a few AU. (The TW 
Hya disc is observed very close to face-on, so the 
line-width cannot be explained by Keplerian rota- 
tion unless the emission arises very close to the star.) 
The profile of the [Neii] line will differ depending 
on whether it arises in a bound disc atmosphere (X- 
ray excitation) or a photoevaporative wind (EUV). 
Unfortunately the observation of this single, face- 
on sour ce does not break the X-ray/EUV degener- 
acy, but lHerczeg et al. (2007) suggest similar obser- 
vations of a sample of discs at different inclination 
angles should resolve this uncertainty. 

Our discussion so far has largely neglected the 
infiucncc of non-ionizing, FUV heating of the disc, 
primarily for reasons of clarity. As we have seen 
in Section 2.2.1, however, FUV heating is capable 
of driving a significant wind from the outer parts 
of the disc, and should be considered in models of 
disc evolution. The much more complex heating 
and cooling problem in this case means that simple 
models in general do not suffice, as the dominant 
heating and cooling mechanisms change with the lo- 
cal disc conditions. The geometrically simpler case 
of FUV heating by e xternal irradiation has been 
modelled in depth by I Johnstone et aL ( 1998h and 
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Adams et alj (|2004l ). and was discussed in Section 
2.2.1. However, to date no self-consistent models 
of FUV ptiotoevaporation by the central star have 
been published. Wo r k in pr ogress was discussed by 
Gorti &: HollenbachI (|2004l ) and lOullemond et al " 



(|2007n . and although the details remain uncertain 
the basic principles seem secure. FUV irradiation 
tends to heat the disc to much lower temperatures 
than EUV heating (typically 100-lOOOK), so FUV 
driven winds tend to be launched from correspond- 
ingly larger radii (> lOOAU). Moreover, the FUV 
heating rate from the central star is likely dominated 



by the accretion-shock emissio n (jMatsuvama et al 



2003 : Alexander et al. I, l2004ah . so time-dependent 
models will behave differently from those of EUV 
photoevaporation. It seems likely that FUV photo- 
evaporation will have a strong influence on the sizes 
of discs around TTs, and may in principle remove a 
significant fraction of the disc mass. However, FUV 
heating is unlikely to affect the disc significantly 
at radii < 50AU, and therefore will not affect the 
EUV-driven clearing processes discussed above. In 
the context of our picture of disc evolution it seems 
likely that FUV-photoevaporation will act through- 
out the evolution, and may shorten the duration of 
the viscous phase of evolution by removing much 
of the disc mass. However, these details remain un- 
certain, and will hopefully be addressed in the near 
future. 



2.5. Summary of gas disc evolution 

In this Section we have discussed the various pro- 
cesses that affect the evolution of gas discs, and 
considered in detail models of disc photoevapora- 
tion. We have discussed external photoevaporation 
in massive clusters, and the application of these 
models to the proplyds observed in the ONC. We 
have shown that models which combine photoevap- 
oration by the central star with viscous evolution of 
the disc predict rapid disc clearing after a long life- 
time, and shown that such models satisfy the obser- 
vational constraints discussed in Section 1. Lastly, 
we have discussed current and future research in this 
area, from both theoretical and observational per- 
spectives. 



3. Dust dynamics and evolution 

We now turn our attention to the solid com 
nent of the disc, generally referred to as "dust 
Dust is thought to make up only around 1% of the 
total mass of protoplanetary discs, but it is crucial 
to the study of disc evolution for several reasons. 
Firstly, the magnetohydrodynamic (MHD) turbu- 
lence that is thought to drive angular momentum 
transport in accretion discs requir es that the disc 
be ionized above some low level ( Gammi^ . 1996 : 
Balbus fc Hawlev . 19981 ). The dominant source of 
ionization in the midplane of protoplanetary discs is 
thought to be ionization of heavy e lements by high- 



energ y photons and cosmic rays (jGlassgold et al 



2000; Fromang et al. I. l2002h . and in this region most 
of the heavy elements are "locked up" in dust grains. 
In addition, free electrons in the disc are readily ab- 
sorbed by dust grains, and consequently an under- 
standing of the distribution of dust in the disc is 
crucial to our understanding of gas disc evolution. 
Secondly, as discussed in Section 1, while dust rep- 
resents only a small fraction of the mass in proto- 
planetary discs, it contributes most of the disc opac- 
ity. As a result most observations of protoplanetary 
discs observe the dust rather than the gas, so ideally 
models which seek to reproduce such observations 
should treat the dust and gas components of the disc 
independently. Additionally, most TT discs are pri- 
marily heated by stellar irradiation, and as the dust 
dominates the disc opacity it is therefore crucial to 
our understanding of the thermal structure of discs. 
Lastly, and most importantly for this chapter, the 
dust in discs is thought to represent the "building 
blocks" of planets. Any theory of planet formation 
must explain the growth of dust from the small, sub- 
lim grains found in the interstellar medium to the 
very large solid bodies and planetary cores we see 
in planetary systems. The evolution and dynamics 
of dust in discs is therefore of crucial importance to 
our understanding of both disc physics and planet 
formation. 

3.1. Basic physics of grain growth 

We first address the basic physics of dust growth, 
using a very simple approach. Much of the discus- 
sion that follows is originally attributed to Victor 



^ In an astrophysical context the term "dust" refers to solid 
material across a very large size range, from sub-/jm sized 
grains up to rocks of metre-size and larger. 
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Safronov, whose 1969 monograph Evolution of the 
protoplanetary cloud and the formation of the Earth 
and the planets stiU forms the basis of most modern 
theories of planet formation. Here I discuss some of 
the simple arguments, first presented by Safronov, 
on the growth of grains in protoplanetary discs. 

We can estimate the rate at which grains grow 
very simply, by considering collisions between dust 
grains in the disc. We assume a typical material den- 
sity of dust particles of pd = Ig cm^"^, and note 
that typical midplane gas densities in TT discs are 
pg ~ 10~®g cm^'^. If the grains are sufficiently small 
that they are very well coupled to the gas by aero- 
dynamic drag (see Section 3.2 below), then the ve- 
locity dispersion, a, of the grains is approximately 
that of the Brownian motion of the particles, which 
can be approximated as 



TOH 



(29) 



If we assume that small particles stick with 100% 
efficiency when they collide, the mean growth rate 
of particle mass is simply given by 



dm 
~dt 



ns^PgZa . 



(30) 



Here s is the radius of the dust particles, which are 
assumed to be spherical, and Z is the dust-to-gas 
ratio (metallicity) in the disc. The mass of a grain 
m = (4/3)7rs'^pd: so the rate of growth of particle 
size can be written as 



ds ,P^r^ 
— 4:!-^Za. 

dt Pd 



(31) 



The typical sound speed in a disc at 1 AU is Cg — 1km 
(which gives a ~ 0.1cm s~^), and if we assume 
a canonical dust-to-gas ratio oi Z — 0.01 we find a 
growth rate of 



ds _A 

— ~ 10 cmyr 

dt 



(32) 



Given that the grains move relative to the gas (see 
below), the velocity dispersion a is likely to be 
significantly larger than we have estimated above. 
(Strictly, this estimate applies only to the small- 
est grains, with s < O.lfim.) Consequently, even if 
we assume that grains do not stick with 100% effi- 
ciency we still expect to see rapid growth of small 
grains, with grains growing to cm-size on a timescale 
< lO^yr. However, once grains reach this size range 
they no longer remain well-coupled to the gas in 
the disc, and this analysis breaks down. Instead we 



must consider the effects of gas-drag forces on the 
motions of the dust grains: forces which produce 
some surprising results. 

3.2. Motion of dust grains 

We now consider the motion of dust particles 
in protoplanetary discs. Such grains are subject to 
gravity, centrifugal forces and aerodynamic drag, 
and so we will now discuss the effects of aerody- 
namic drag on dust grains in discs. A spherical 
particle of radius s, moving at a velocity v relative 
to gas of density pg experiences an aerodynamic 
drag force which opposes its motion: 



1 



(33) 



This expression has three terms: the cross-sectional 
area of the grain, tts^; the ram pressure exerted on 
the grain, Pgv'^; and the drag coefhcient Cd. The 
drag coefficient depends on the velocity of the grain 
relative to the gas, and also on the size of the grain 
relative to the mean-free-path o f gas molecules A 
(IWhippld . I1972I : IWeidenschillh^ . Il977l) . For small 
particles (s < 1cm), where the size of the particles is 
less than A [formally where s < (9/4)A] we are in the 
Epstein regime, and the drag coefficient is given by 



8 vth 
3 V ' 



(34) 



where Uth = (S/tt) ' Cg is the mean thermal velocity 
of the gas. For larger particles we instead find our- 
selves in the Stokes' regime, where the drag coeffi- 
cient depends on the Reynolds number of the flow. 
In this regime Cd is typically approximated by a 
piecewise fun ction of the Reynolds n umber, such as 
that given in IWeidenschillingl ( 1977 ) . 

It is useful to define the stopping timescale (or 
drag timescale), tg, which is the timescale on which 
frictional drag will cause an order-of-unity change 
in the momentum of the dust grain (rela tive to the 
gas) . Here we define the stopping time a43° 



(35) 



so in the Epstein regime the stopping time is there- 
fore given by 



^''Note that different authors tend to define the stopping 
timescale sUghtly differently, often with different numerical 
pre-factors. Consequently care must be taken when compar- 
ing different works, as the exact definitions of tg may change. 
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Pd s 

Pg Vth 



(36) 



If we adopt values typical of conditions in a proto- 
planetary disc (e.g., pg = 10~^g cm~^, vth — Ikni 
s~-^) and the standard dust density of pd = Ig cni~^, 
we see that the stopping time for /xm-sized grains 
is of order seconds. This is much smaller than the 
typical dynamical (orbital) timescale, so we see that 
small grains are extremely well-coupled to the gas. 
However, metre-size particles have ts ^ and are 
only marginally coupled to the gas, and much larger 
bodies are mostly unaffected by gas drag. 

3.2.1. Settling 

The simplest motion of dust grains we can con- 
sider is that of vertical settling in an isothermal gas 
disc. We estimate the settling timescale by equating 
the opposing forces of drag and gravity. For vertical 
displacements z <^ R, the vertical component of the 
gravitational force on a grain is simply 



(37) 



where VIk = {GM^/ R^Y/'^ is the Keplerian orbital 
frequency. If we equate this with the drag force (in 
the Epstein regime) and re-arrange we find a typical 
settling velocity of 



wth Pg,{z) 

we can then estimate the settling timescale as 
, _ z _ Vt\, pgjz) 1 

f-scttlc — — ^9 

^settle "k Pd S 



(38) 



(39) 



For /im-sized particles in a typical protoplanetary 
disc the settling timescale is ~ lO^yr, but larger 
grains are expected to sediment towards the mid- 
plane in much shorter timescales. Note also the de- 
pendence on the gas density: in a vertically isother- 
mal disc Pg oc exp(— z^/2i/^), so we expect parti- 
cles to sediment out of the upper layers of the disc 
rapidly, but settle more slowly at smaller z. 

In principle vertical settling will increase the mid- 
plane dust-to-gas ratio, and can result in increased 
dust collision (and growth) rates. However, this 
simple analysis considers only a static, laminar ver- 
tical disc structure, and real TT discs are thought 
to be turbulent. Settling of dust in a turbulent disc 
is a much more complex problem, as the turbulence 
can lift the grains to high z on the eddy timescale 
of the turbulence, which is much shorter than the 
expected dust settling timescale. Several recent 



studies have used numerical simulations to study 
the effects of MHD t urbul e nce on dust settlin g (e.g. . 



Johansen fc Klahii 120051: ICarballido et all . l2005t 



Turner et al.l. 2009), and while some details remain 



contested, all the models agree that turbulence does 
result in much slower settling than predicted above. 
Indeed, it seems likely that the vertical distribution 
of small grains in discs is sustained by some sort of 
equilibrium between turbulent diffusion and sedi- 
mentation. (For further discussion of the effects of 
turbulence see the chapter by Klahr in this volume.) 

3.2.2. Radial drift 

The equations governing the radial drift of 
dust in protopl a netary disc were first derived by 
Weidenschillind ( 1977); here we follow the more 
recent analysis of iTakeuchi &: LinI (|2002ll2005f) . We 
first derive the relationship between the rotational 
velocity of the gas, v^^g and the Keplerian orbital 
velocity vk = {GM^/ R)^/"^ . We consider the radial 
component of the equation of motion thus: 



R 



GM, 

i?2 



1 dp 

'p^dR 



(40) 



Here the left-hand side is the centrifugal accelera- 
tion, which balances the accelerations due to grav- 
ity and gas pressure on the right-hard side. If we as- 
sume that the gas pressure at the disc midplane can 
be approximated as a power-law 

P (X i?-" (41) 

and assume a locally isothermal equation of state 



p = pci 



(42) 



then the pressure gradient term in Equation 40 be- 
comes 

(43) 



1 dP _ c; 
Pg dR R 



If we substitute this into Equation 40 and multiply 
by R, we find that 



where the term 



V 



(44) 



(45) 



K 



denotes how sub-Keplerian the gas is. Typical 
power-law discs have values of rt ~ 2.75-3, depend- 
ing on the choice of viscosity and temperature laws. 
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Therefore, if we assume that H/R = c^/vk = 0.05 
and n = 11/4 (consistent with an a-disc with 
T (X we see that the gas in a protoplanetary 

disc is typically sub-Keplerian by 



(1 - r;)i/2 ^ 0.0034. 



(46) 



This may seem like a small fraction, but bear in 
mind that the Kepler speed for a IM© star at lAU 
is ~ 30km s~^. Consequently, we expect particles 
orbiting at the Kepler velocity to be subject to a 
very strong headwind, ~ 100m s~^. 

We now consider the radial drift of a dust particle 
of radius s orbiting in such a gas disc. Before we do, 
however, it is useful to define a dimensionless stop- 
ping time Tg. This is simply the stopping timescale 
ts divide d by the orbital timescale, and here we de- 
fine Ts aJ^ 



Ts — isf^K — t. 



Vk 

R 



(47) 



If we now consider the radial and azimuthal equa- 
tions of motion (EoMs) for a dust grain, we see that 



r 



-(Wr,d 



dt 



t 



(48) 



(49) 



Here the subscripts r and ^ denote the radial and 
azimuthal components of velocity respectively, with 
the additional subscripts g and d used to distinguish 
the gas and dust velocities. In the radial EoM, the 
first term is the centrifugal acceleration, the second 
term the acceleration due to gravity, and the third 
term the frictional drag force (which opposes the 
motion of the dust grain, and is zero if the dust and 
gas move at the same velocity). In the azimuthal di- 
rection the only acceleration is that due to the drag 
force. We then make the simplifying assumption that 
the dust particles spiral in on approximately circu- 
lar orbits (i.e., the specific angular momentum of a 
particle is always close to Keplerian), so that to first 
order 



(50) 



Consequently we can simplify the azimuthal EoM 
by noting that 

■^{rv^^d) ^ v^^d-^irvK) = ^«0,dWK , (51) 



Note again that exact definitions of Ts differ in the liter- 
ature, commonly by factors of 2tt. 



and by substituting into Equation 49 and re- 
arranging we find that 



ItsVK 

2 r 



1. 



:TsVrd 



(52) 



We then make the second simplifying assumption 
that, to 0{{H/Rf), there is no net acceleration in 
the radial direction (i.e. dv^^d/dt ~ 0). If we re- write 
the Keplerian velocity in terms of the gas orbital 
velocity and rj thus 



n^R = 



(53) 



R R ' 

we can substitute for the appropriate term in the 
radial EoM to find 



-'K 



(Wr,d - Wr,g) = . 



(54) 



We expand the first two terms as a difference of 
squares (noting that w^^d + v^.g ~ 2wk), and substi- 
tute the difference between the azimuthal velocities 
of the dust and gas from Equation 52. If we then 
multiply by vk/R and re- write in terms of Tg, we 
find that 



(55) 



and therefore the drift velocity of the dust is given 
by 



Wr,d 



rjVK 



(56) 



The form of this expression is shown in Fig. 7. The 
radial velocity of the gas (due to viscosity) is typi- 
cally very small, and if we neglect this motion we see 
that Equation 56 is symmetric in Tg. The limit Tg <C 
1 corresponds to very small particles (s < 1/^m), and 
in this limit the grains are so well-coupled to the gas 
that they simply migrate inwards at the same speed 
as the gas. By contrast, in the limit Tg ^ 1, which 
corresponds to very large particles (s > 1km), the 
particles have so much inertia that they are hardly 
affected by gas drag at all, and simply orbit at the 
Kepler speeqjfj. The dynamics of bodies in this size 
regime, commonly referred to as "planetesimals" , 
are in fact dominated by gravitational interactions 
between the planetesimals, and it is thought that 
gravitational focusing increases the collision rates 



Note, however, that for Ts > 10 we are in the (non-linear) 
Stokes regime, rather than the Epstein regime, so the simple 
linear relationship between particle size and stopping time 
(Equation 36) no longer holds. 
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Drift velocity in flaring disc at 1AU 



are not seen (e.g.. [Andrews fc Williama . 2005) but 




Fig. 7. Inward drift velocity, in units of the Kepler speed, for 
a typical TT disc at lAU (see Equation 56). The departure 
from symmetry is due to the inward gas velocity, which is 
computed for different ct parameters at an accretion rate of 
10~*MQyr~^. Note the very large drift velocity for grains 
with Ts ~ 1: for a IMq star, IQ-^dk 3000cm s-1. 

between such bodies sufficiently th at they ca n grow 
to planetary sizes very rapidly (e.g.. lSafronovlll96S ) . 

However, we also see that from Fig. 7 that at the 
peak of the curve the inward drift is very rapid, 
reaching > lOOOcm s~^ in a disc around IMq star 
at lAU. At such velocities dust particles will spi- 
ral in to the central star in ~ lOOyr, and for typi- 
cal TT disc parameters, Tg = 1 corresponds to sizes 
of s ~ 50-lOOcm. We see, therefore, that in order 
to form planetesimals, particles must grow through 
the range of sizes that correspond to Tg 1 on very 
short timescales, and the conclusion is that planetes- 
imal formation must occur very rapidly is unavoid- 
able. This is commonly referred to as the "metre- 
size barrier" , and much modern research on p lane t 
formation focuses on how to solve this probleir P^ . 

In recent years several models have begun to 
combine models of gas evolution (similar to those 
discussed in Section 2) with radial drift of dust, 
with a vie w to modelling dust e volution over Myr 
timescales. Takeuchi et al.l (l200d) combined the two- 
fluid model ofl Takeuchi fc LinI (|2002l . l2005l) with the 



photo evaporation/ viscous model of IClarke et al 



(|200lh ;o study the dynamics of mm grains, and 
found that, if unreplenished, the mm grain popu- 
lation in TT discs is accreted on to the star before 
the gas disc is cleared. Dust-poor discs of this type 



Cynically one could say that if we did not know of the 
existence of planets, then the metre-size barrier would be a 
very robust argument as to why planets cannot form via the 
coUisional growth of smaller particles! 



mm- to cm-size pa r ticles are observed ([Testi et al 



Takeuchi et al 



2003t IWilner et al.l. l2005t iRodman et all . l2006f) . 

(|2nn5h 



so 

concluded that some addi- 
tional effect, such as grain growth or a reduced 
radial drift veocity, must slow the decline in t he mrn 
grain population. [Alexander fc Armitagd (|2007l ) 
extended the study of Takeuchi et al.l ( 2005 ) to in- 
corporate more realistic photoevaporation models, 
and also considered a range of grain sizes. However, 
without considering grain growth their model was 
restricted to timescales of < lO^yr, so they focused 
on the behaviour of the dust during the clearing 
phase of the ev o lution (see Section 4). More re- 
cently, iGaraudI (|2n07h considered models which 
include grain growth and dynamics, and found that 
the long-term evolution of the dust population de- 
pends very sensitively on the sizes of the grains 
that dominate the mass reservoir at lar ge radius (a 
result also found bv iBrauer et al.l . l200 /). Models of 
this type are becoming increasingly sophisticated, 
and provide a promising avenue for future research. 



3.3. Dust growth and planetesimal formation 

We have seen in Section 3.1 above that the growth 
of small dust grains can, in principle, occur very 
rapidly. More realistic models than those presented 
above take into account more realistic physics of 
collisions and sticking, and in recent years much 
progress has been made, through both theoreti- 
cal calculations an d labo r atory experiments (e.g. , 
Blum fc Schrapleil 12004 iDullemond fc Dominikl . 
20051, see also the review bv iDominik et al.ll2007t ). 
Sticking efficiencies have been found to depend on 
a large number of factors, such as the shapes of 
the colliding bodies, the porosity of the grains, or 
the velocity of the collision. In general it is found 
that sticking efficiencies tend to be high (10-100%) 
at small grain sizes, and indeed the growth of 
small grains is so rapid that many models of dust 
growth require some form of growth-destruction 
equilibrium in order to sustain the observed popu- 
lation of /.tm-sized grains over Myr timescales (e.g 



Dullemond fc Dominik . l2005h . However, sticking 



becomes much less efficient once grains reach the cm 
size range, resulting in much slower growth at these 
sizes. Moreover, in this size range grain collisions 
can become destructive at speeds > 10m s~^, so the 
dynamics of the gas-grain interact ion become im- 
portant for growth processes also (jPominik et al 
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20071) 



Many popular mechanisms for the formation of 
planetesimals are rooted in the so-ca. l led G oldreich- 
Ward instabihty (jColdreich fc Wardl . [l973h . In this 
scenario the vertical settling of small dust grains 
steadily increases the dust-to-gas ratio at the mid- 
plane, until eventually the thin dust sub-disc be- 
comes gravitationally unstable and fragments. We 
can estimate the conditio ns at which this occurs by 
a simple application of the Toomre ( 19641 ) criterion, 
setting 



Q = 



ttGE 



= 1. 



dust 



(57) 



where a is the velocity dispersion in the dust layer 
and Sdust is the surface density of the dust layer. If 
we assume that Sdust = Sgas/100 ~ lOg cm~^ and 
consider a IM© star, we see that a velocity disper- 
sion of (7 < 10cm s^^ is required in order for gravita- 
tional instability to fragment the dust layer. Given 
that the gas sound speed is typically ^ 1km s~^, we 
see that the dust sub-layer must be very thin indeed 
in order to become unstable. However, this sugges- 
tion is attractive because it allows for the formation 
of planetesimals directly from small (s < 1mm) par- 
ticles, and thus bypasses the size ranges what are 
most susceptible to radial drift. 

Unfortunately we can easily show that, in its 
simplest form, the Goldreich-Ward mechanism does 
not result in planetesimal formation. Wc have seen 
above that the instability requires that the dust 
layer be approximately 10** times thinner than the 
gas disc, and consequently at the disc midplane the 
local dust density dramatically exceeds the local 
gas density (by a factor ^ 100). In this layer we 
therefore expect the solid material to dominate the 
dynamics, so gas pressure forces are negligible and 
the dust and gas both orbit at the Kepler speed. 
Above this layer, however, the gas is significantly 
sub-Keplerian due to gas pressure (as discussed 
in Section 3.2.2 above), and there is therefore a 
large velocity shear in the vertical direction. This 
shear is Kelvin-Hclmholz unstable, and the result- 
ing turbulence prevents cr fro m becoming small 
enough for instability to set in ( Cuzzi et al.l . [l993l. 
More recently, several sets of authors have argued 
that planetesimal formation can occur in this man- 
ner, if the local dust-to-gas ratio is significantly 
enhanced from the c a nonical value o f 1:100 (e.g., 
Garaud fc Li^ . l2004l) . lYoudin fc ShT] (l2002l) sug- 
gested that radial drift can result in large local 
enhancements in the dust-to-gas ratio, leading to 



instabihty, while IXhroop fc Ballvl (|2005l ) suggested 
that photoevaporation by external irradiation 
(which preferentially removes gas from the disc) can 
likewise lead to grav itational instability i r i the d ust 
laye r. More recently, Johanscn fc Youdin J 20071 see 
alsolYoudin fc Goodman, 2005; Youdin fc JohansenI 
2007ir showed that radial and azimuthal drift of 



particles in a gas disc can lead to the so-called 
streaming instability, and that the resulting struc- 
ture in the gas disc can significantly increase the 
concentration of solid bodies. They find that par- 
ticles which are marginally coupled to the gas are 
most prone to clumping in this manner, and also 
find that the net radial drift of these particles is 
significantly reduced. This mechanism provides an 
attractive solution to the metre-size problem, and 
represents one of the mo st promising me ans of 



2003). 



forming planetesimals ( Johansen et al. 

In a similar vein, much recent research has in- 
voked structure in the gas disc in order to accelerate 
planetesimal formation by increasing collision rates. 
A variety of similar models have been proposed, but 
all make use of the same underlying physics. We saw 
in Section 3.2.2 that gas pressure gradients in pro- 
toplanetary discs result in radial migration of dust 
grains, but our treatment of the problem assumed 
a simple disc where the uniform pressure gradient 
invariably makes the gas slightly sub-Keplerian. 
This results in inward migration of dust grains, but 
in general (as can be seen from Equations 40 fc 56), 
the dust will in fact mov e radially t oward s any local 
pressure maximum (;e.g.. lWhiDr)y . ll972h . If the gas 
disc has some local structure, such as spiral density 
waves, the gas orbits at slightly super-Keplerian 
speeds at the inner edge of the pressure enhance- 
ment. The dust in this region instead migrates out- 
wards, towards local pressure maxima, which can 
result in significant enhancements in the local dust- 
to-gas ratio. Moreover, the particles which are most 
susceptible to this behaviour are the same particles 
that suffer the most severe radial drift, so models of 
this type offer an attractive solution to the metre- 
size problem. This behaviour has been noted by 
a number of authors, an d several different scenar- 
ios have been proposed. Rice et al] (2004, 2006a ) 
demonstrated that the spiral density waves induced 
by gravitational instability in the gas disc result 
in strong local enhancements of the dust-to-gas 
ratio in metre-size particles, and suggest that the 
enhancement may be large enough to allow plan- 
etesimal formation by gravitational interactions 
between the solid bodies. iPaardekooper fc Mellema 
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observed similar behaviour in the spi- 
ral waves produced by the presence of a planet in 
the disc, and suggested that the presence of one 
planet may act to accelerate t he formation of sub - 
seq uent planets. In additi o n. iRice et alJ (l2006bri 
and Alexander &: Armitag^ ( 2007L see also GaraudI 



2007t ) both suggest that local pressure maxima can 
provide a "size sorting" mechanism in discs with 
gaps or holes, as grains with Tg — 1 tend to clump 
in local pressure maxima, while smaller particles 
remain coupled to the the gas flow and larger par- 
ticles orbit independently of the gas. All of these 
models suggest that particle collision rates can be 
significantly enhanced if the gas disc has sufficient 
structure, but whether or not this occurs in real 
systems remains unclear. In real discs any, all, or 
some combinations of these various mechanisms 
may apply, but the details of planetesimal formation 
remain rather poorly understood. 



3.4. Summary of dust evolution 

In this Section we have discussed how dust grains 
grow from small, sub-^m sizes to cm-size and larger. 
We have considered the effects of aerodynamic drag 
forces on dust particles, and derived expressions for 
both the vertical settling and radial drift that result. 
Particles of lO-IOOOcm are especially susceptible to 
radial drift, and we have seen that this "metre-size 
barrier" provides a significant challenge for theories 
of planetesimal formation. We have discussed vari- 
ous mechanisms for forming planetesimals, but have 
noted the problems associated with all of them. How 
to form planetestimals remains the biggest challenge 
for modern research in planet formation, and many 
of the details of this problem still elude us. 



4. Transitional discs 

In recent years an increasing literature has built 
up on so-called "transitional" discs. This small but 
important class of objects have properties between 
those typical of CTTs and WTTs, and as such are 
thought to represent discs which are observed during 
the (short) clearing phase of their evolution. These 
objects are crucial to our understanding of disc clear- 
ing but, as the clearing process is rapid (see Section 
1.2), the objects are rather rare. To date around 20 
such objects are known, but this sample is expand- 



ing rapidl}{^. In general it is believed that these ob- 
jects show characteristics indicative of some degree 
of inner disc clearing, and different authors have in- 
voked a variety of different physical mechanisms to 
explain these "inner holes" . In this section I attempt 
to summarize this rapidly-evolving subject, and try 
to highlight some key areas for future research in 
both observational and theoretical fields. 

4.1. What is a "transitional" disc? 

Despite appearances, this is not a trivial question. 
As mentioned above, a rough definition would en- 
compass all objects with properties between those 
typical of CTTs and WTTs but, given the broad 
range of observations used to study protoplanetary 
discs, it should come as little surprise that a ro- 
bust definition has so far been difficult to establish. 
In general, objects have been classified as "transi- 
tional" based on their SEDs, and a loose definition 
would be that such objects typically show a fiux 
deficit in the near- and mid-infrared, when com- 
pared to "ordi nary" CTTs . The first such objects 
were found by IStrom et all (|l989f ). who coined the 
term "discs in transition" and suggested that the 
observed SEDs were indicative of some degree of in- 
ner disc clearing. Around 20 such objects are now 
known, but the advent of Spitzer is leading to a rapid 
increase in this sample size. However, there is little 
consistency in the definitions used to classify these 
objects, in terms of either nomenclature or observa- 
tional characteristics. Recently, objects with similar 
prop erties have va. r iously been called "transi tional" 
fe.g-. ICalvet et all. l2005l: iNaiita et al.l. l2007l). "pas- 
sive" (iMcCabe et al.l.l2006l) . "a naemic'' (iLada etaT 



2006[ ) or "cold" (jBrown et all 120071 ). and a similar 
variety of observational selection criteria exist in the 
literature. 

The best-studied transition objec ts are th e four 
nearby sourc es TW Hya (e.g., C alvet et al.l . [2002; 
Wilner et all |2005| . found in the TW Hya acco- 
ciation at a distance of ~ 60pc), CM Aur (e.g.. 



Rice et al.l. 120031: ICalvet et al 



fe.g.. ICalvet et al. 



PC ), uiy i 

1 limi, 



.T...*^,— , DM Tau 

- I20n5^ and CoKu Ta u/4 (e.g., 

Forrest et al.l . l2004Hd'Alessio et al.l . l2005al . all found 
in the Taurus- Auriga cloud at a distance of ~ 
140pc), and in essence these four objects have come 
to define the modern meaning of "transitional disc" . 
All four show little or no excess emission above 



Around a dozen new papers on transitional discs appeared 
while this article was being written! 
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their stellar photospheres in the near-infrared, but 
show excesses comparable to conventional CTTs 
at wavelengths longer than ~10-20/j,m. (A typical 
spectrum, that of CoKu Tau/4, is shown in Fig. 8). 
Detailed modelling has shown that the SEDs of 
these objects are best-fit by discs with large inner 
holes: that is, discs that are (mostly) cleared of dust 
at small radii, <5-10AU, but appear to be "nor- 



mal" CTT disc s at larg e r radii (ICalvet et al l . 12002 



d'Alessio et aP . l2005at ICalvet et all . l2005[ ). The 



best-fitting hole radii are typically a few to a few 
tens of AU, and do not appe ar to sca le simply with 
any known stellar propertie d ^^H ^^1 . A large (and 
increasing) number of different mechanisms have 
been proposed to explain these inner holes, with 
different studies variously invoking the presence 
of a planet (e.g., Rice et all . 12003 : Quillen et al 



2004) photoevaporation (e. g.. [Alexander et al 



2006bD . dust e volution (e.g. IWilner et all . 120051 ). 
photophoresis (jKrauss et al 1. l2007h or enhance- 
ment of the magnetor o tation al instability (MRI, 
Chiang &: Murray-ClavL 120071 ) to explain some or 
all of the observed properties. 

However, when we look beyond the SEDs we see 
that observations of transitional discs show remark- 
able diversity. For example, the observed (stellar) 
accretion rates and disc masses for these four objects 
alone differ by more than 3 orders of magnitude. 
Moreover, the structures of these four discs seem to 
differ significantly within the holes. GM Aur and 
DM Tau are both CTTs with significant gas accre- 
tion (^ lO^^M0yr^^) observed at the stellar sur- 
face, but while GM Aur shows a weak near-infrared 
excess, suggestive of some optically thin dust close to 
the star, DM Tau appears to have a much "cleaner" 
hole in its dust disc (Calvet et al.. 2005). TW Hya is 
accreting at a much lo wer level (~ 4x lO~^°M0yr~^, 
Muzerolle et~al]l2000f ). but interferometric observa- 
tions of the inner disc suggest the presence of a pop- 



Note, however, that the sizes of the holes are not espe- 
cially well-known. Recent mid-IR inte rferometric observa- 
tions of TW Hya l|Ratzka et al.l. l2007t ) measured a signifi- 
cantly smaller inner hole radius than had pr eviously been 
measured from the SEP llCalvet et al.L l2002h or from mm 
interferometry l lHuehea et al.l. 120071) 7 and the origin of this 
discrepancy is not yet understood. 

1^ Recently, IChiang fc Murrav-Clavl l|2007l ') have claimed 
that the observed hole sizes for accreting transition discs 
correlate with the measured accretion rates, with M oc 
arj^^j^/M* . The small number of sources (three) and large 
error bars make the significance of this correlation difficult 
to assess, but future observations should be able to either 
confirm or reject this hypothesis with high confidence. 



ulation of small dust grains that must be continually 
replenis hed, presumably by the destruction of larger 
bodies ( Eisner et al.l . l200d) . CoKu Tau/4, by con- 
trast, is a WTT with no detectable accretion, and 
shows no evidence for significant amounts of gas or 
dust within its hole. 

As discussed in Section 1, protoplanetary discs 
are generally extremely optically thick in the near- 
and mid-infrared, so the bulk of the disc emission at 
these wavelengths is ex pected to arise in a thin laye r 
on the disc surface (e.g., Chiang fc Goldreiclil . ll997 ). 
To first order, therefore, the infrared BED only tells 
us about the presence or absence of dust at a given 
temperature in the disc. All of the models for disc 
clearing mentioned above result in similar infrared 
SEDs, and as such infrared observations alone are 
of limited use in distinguishing between models. In- 
deed, no single object is unambiguously associated 
with any particular model. Consequently it seems 
likely that a statistical approach, selecting objects 
on the basis of their SEDs and then investigating 
other disc properties, may be the most promising 
way of learning more about these objects. 

The only such st udy published to date is that of 
Najita et al.l ( 2007t l. who used the Spitzer I n frared 



Spectrograph (IRS) survey of ^ Furlan et al. ( 2006h 
to identify a sample of 12 transitional discs in the 
Taurus- Auriga cloud. Taurus is uniquely suited to 
a statistical study of disc evolution, as the stellar 
properties of all of the sources are well known and 
the SCUBA survey of I Andrews fc Williams] (|2005h 
provides a uniform sample of disc masses for al- 

lloo^) 



Naiita et al. 



most all of the known TTs 
utilised a demographic approach, and showed that 
their transitional objects show systematically lower 
accretion rates than CTTs of the same disc mass. 
This strongly suggests that objects with transi- 
tional SEDs do represent an evolved state of CTTs, 
in agreement with previous studies of individual 



objec ts (e.g., Id'Alessio et al.l . l2005al: ICalvet et al 



2005). However, the work of Nanta et al.l 



et al-L 

mm 



introduces a potential selection bias by selecting 
objects based on their infrared SEDs alone. This 
is not necessarily problematic but, as the authors 
discuss, is not an ideal method of selecting a sample 
to compare with theoretical models. Consequently, 
I will now attempt to outline a theoretically moti- 
vated definition of what a transitional disc is. 

All of the theoretical models mentioned above 
(with the exception of dust evolution) produce deep 
holes in the inner (dust) disc, and a sharp transition 
between the inner, cleared region (the hole) and the 
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outer disc. Consequently I propose that the selec- 
tion criterion used to define a transitional disc for 
the purposes of comparison with theoretical models 
should be one that unambiguously requires a hole 
in the dust disc. Consequently objects should be 
classified as "transitional" based on the following 
criterion: 

Transition objects have discs which show optically 
thin emission at shorter wavelengths AND optically 
thick emission at longer wavelengths. 

Objects which satisfy this criterion must have ra- 
dial disc structures which are clearly inconsistent 
with those typical of CTTs, and essentially must 
show strong evidence for a hole or gap in the inner 
regions of the disc. In fact, for reasons of clarity, it 
may well make sense to class such objects as "in- 
ner hole" discs, essentially defining a sub-set of the 
more loosely identified transitional discs. 

An example of the application of this criterion is 
shown in Fig. 8. Shown in the figu re is the median 



SEP for CTTs in Taur us- Auriga (jd'Alessio et al 



1999l:lFurlan et al.l . l2006h . as well as the SEDsoftwo 



transi tion disc taken from the sa mple of lNajita et al 



(|2007l ): CoKu Tau/4 and FQ TarLj. Also shown is 
the spectrum of a stellar black-body (with a tem- 
perature of 3200K, typical of an early M star), and 
that of a simple disc model. The disc model assumes 
power-law surface density, and temperature, 

T(i?), profiles as follows: 



Y.{R) = 100 



T{R) = 300 



R 



(lAU 



R 
lAU 



gcm~2 ^ i? < lOOAU , (58) 



-1/2 



K , R< lOOAU. (59) 



The minimum disc temperature is taken to be lOK, 
and the inner edge of the disc (i.e., the dust subli- 
mation radius) is fixed at the point where T{R) ~ 
1500K. The outer disc radius of lOOAU is adopted 
merely for convenience, and does not affect the emis- 
sion in the infrared. The disc SED is then computed 
as 




Fig. 8. Infrared spectra of transition discs. The median SED 
of Taurus CTTs is shown in blu e. The photometric median 
SED from d'Ale ssio et al. I lll999h is shown as p oints with er- 
ror ba rs, and the Spitzer IRS median SED from lFurlan et aP 
is shown as solid and dashed lines: plotted in each 
case are the median and the quartiles of the observed distri- 
bution of SEDs. The solid black lines show the stellar black 
body and simple disc model described in the text, and the 
dashed and dotted lines respectively show the optically thick 
and optically thin (r = 0.05) limits of the disc model. The 
SEDs of CoKu Tau/4 and F Q Tau are plotted in green and 
red respectively (data from iFurlan et al.l l2006l . kindly pro- 
vided in electronic form by Joan Najita): the lines are IRS 
spectra; the diamonds and triangles denote IRAS detections 
and upper limits respectively. 

where d is the Earth-star distance (taken to be 
140pc), B^{T} is the Planck function and i = 60° is 
the inclination angle of the disc. The optical depth 
is evaluated as 



= Ku^{R) , 



(61) 



ass uming a simple power-l aw form for the dust opac- 



ity (jBeckwith et al.l . ll990f) 
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(62) 
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'^'^ These particular objects are chosen simply for illustrative 
purposes, as they have similar spectral types (Ml. 5 and M2 
respectively) and luminosities. FQ Tau is also known to be 
a binary system, but this does not affect this illustrative 
discussion. 



and the final SED is then computed by adding the 
contribution from the disc to the stellar black-body. 

Also shown in Fig. 8 are the optically thick and 
optically thin (r = 0.05) limits of the disc model. 
In order to satisfy the criterion above, objects must 
have SEDs comparable to or below the optically thin 
limit at short wavelengths, but consistent with the 
^gQoptically thick limit at longer wavelengths. [In prac- 
tice, one could use the (normalised) quartiles of the 
median SED and the 8~10/im stellar (photospheric) 
flux to represent the optically thick and thin lim- 
its respectively.] We see from Fig. 8 that the spec- 
trum of CoKu Tau / 4 shows essentially no excess over 
the stellar photosphere at wavelengths shorter than 
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8/ini, but is consistent with an opticaUy thick disc at 
wavelengths longer than ~ 17/im. By contrast, the 
IRS spectrum of FQ Tau shows a weak excess across 
the entire IRS band, and is consistent with an ap- 
proximately constant optical depth of r ~ 0.1. The 
IRS spectrum of FQ Tau, therefore, does not require 
a genuine hole to be present in the disc, but could in- 
stead be fit by a disc with lower-than-average optical 
depth (possibly due to dust settling and/or growth). 
In this case it seems likely that the disc does emit in 
the near-infrared, but that the relatively weak disc 
emission at such wavelengths results in a negligible 
excess over the stellar photosphere. 

However, the IRAS 60/xm observation of FQ Tau 
does show a significant excess, an order of magni- 
tude above the optically thin limit, and consistent 
with the simple disc model at this wavelength. This 
highlights what is potentially an important selec- 
tion effect: that the wavelength range covered by 
the SpitzerlYlS is only sensiti ve to fairly small holes, 
up to a few tens of AU in siz Any objects with 
larger inner holes (> 50AU in size, although the cut- 
off depends strongly on stellar spectral type) will be 
defined as Class III objects on the basis of their IRS 
spectra alone, as the disc will only emit at wave- 
lengths longer than those observed with Spitzer. 
This highlights the need for multi-wavelength ob- 
servations when identifying transitional discs, as 
techniques which are only sensitive to inner holes in 
a particular size range may bias subsequent inter- 
pretation of the data. Consequently, I propose that 
future statistical studies identify inner hole sources 
using data that covers the widest possible range of 
wavelengths, as only in this manner will unbiased 
samples be guaranteed. 



in th e disc (e.g.. Rice et al. . 2003t Quillen et al 



2004[ ). and gap-opening by the combined actions 



of photoevaporation and viscosity (as discussed in 
Section 2). I restrict myself to considering these 
two models first for reasons of length, but also 
because these are the only models which tell the 
"co mplete story" . Mech anisms such as photophore- 
sis ( Krauss et al. LJ2OO?) or M RI enh ancement by 
X-rays (jChiang fc Murrav-ClavL 120071 ) may well ex- 
plain some of the observed properties of discs with 
inner holes, but these models require the action of 
some other mechanism to create the inner hole in 
the first place. 



4.2.1. Gap- opening by planets 

The planet-disc interaction is discussed in much 
greater detail in the chapter by Klahr (this volume), 
but here I quickly review the criteria for a massive 
planet to open a gap in a disc. If a planet is suffi- 
ciently massive then the tidal torques it exerts on 
the disc can op en a gap close to the planet (e.g., 
i Takeuchi et al.l . :1996). and models of discs with such 
gaps predict SEDs consistent wit h those of several 
observed tr ansition objects (e.g.. iRice et al.l 12003 : 
Quillen et a l. , 2004) . Two criteria must be satisfied 



for a planet to open a gap in a disc: firstly, the tidal 
(Lindbland) torques exerted by the planet on the 
disc must overcome the local viscous torques in the 
disc; and secondly, the Hill radius (i.e., the region 
where the gravity of the planet is greater than that 
of the star) must be comparable to the disc thick- 
ness. In the case of protoplanetary discs these crite- 
ria are approximately equivalent, and in an a-disc 
the app roximat e crit erion for gap-opening is (e.g.. 



Takeuchi et al..ll99i 



4.2. Distinguishing between models of disc clearing 

As mentioned above, a large number of models 
have been proposed to explain some or all of the 
observed properties of various transitional discs, 
and here I attempt to distinguish between some 
of these models. If we use the selection criterion 
discussed above to select objects with inner holes, 
then it seems likely that the properties of these 
discs should reflect the manner in which the holes 
were formed. I will focus on two mechanisms in par- 
ticular: gap-opening by a massive planet embedded 



> 



Given the large IRAS beam size this large excess 
could be due to a companion, but again this is not important 
for this illustrative discussion. 



vl/2 



(63) 



Here q is the ratio of the planet mass to the stellar 
mass, and the subscript p denotes properties of the 
planet. For typical parameters we find that planets 
with masses > 0.5Mjup will open a gap in the disc. 
(For a more d etailed s tudy of the gap-opening cri- 
terion see Edgar et al.l . [2007a ') . Moreover, the crite- 
rion for gap-opening is effectively independent of the 
disc surface density (and therefore the disc mass), 
so a sufficiently massive planet will be able to open 
a gap in a disc of almost any mass. 

Once the planet has opened a gap in the disc the 
tidal and viscous torques will clear the inner part 
of the disc, creating an inner hole (jQuillen et al 
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2004 IVarniere et all . l2006l ). Over long timescales 
the planet will then undergo Type II migra- 
tion, but if the planet is sufficiently massive this 
times cale can be compara ble to the disc lifetime 



Properties of inner hole sources 



(e.g., ISver fc Clark3 . Il995h . However, the planet 
forms somewhat of a "leaky barrier" , so unless the 
planet is extremely massive (> lOMjup) we expect 
some accretion of material, both on to the planet 
and beyond (e.g.. iLubow et ah . 19991 ). The details 
of the planetary accretion process in discs with 
gaps rema in somewhat uncertain (se e, e.g., the dis- 
cussion m iLubow fc d'Angelol . l2006h . but it seems 
likely that most planet-induced gaps will have some 
accretion across the planetary orbit, and therefore 
some accretion on to the central star. 



4.2.2. Planets or photoevaporation? 

We have seen above that planet-induced gaps 
tend to result in continued accretion interior to the 
planet, and can occur in discs with a wide range 
of masses. By contrast, we saw in Section 2 that 
photoevaporation opens a gap in the disc only 
once the accretion rate has dropped to a low level 
(that of the photoevaporative wind). If we then 
assume some viscosity law in the disc this sets an 
upper limit on the masses of di sc with "photoe- 



vapora t ed" holes (~ SMjup, e. g . [Alexander et al 



2006bt [Alexander fc Armitagj . |2007|) . Moreover 



once viscosity has drained the inner disc we ex- 
pect essentially no accretion on to the central star, 
as the direct wind clears the disc on a timescale 
shorter than the local viscous timescale (see Section 



2.3.3) . In li ght of these fa c ts. [ Alexander fc: Armitage 
(|2007t ) and iNaiita et all (|2007l ) independently pro- 
posed that surveys of disc masses and accretion 
rates can be used to distinguish between these two 
mechanisms for inner disc clearing. 

This is dem onstrated in Fig. 9, wh i ch com bines 
the model s of Alexander fc Armitagd (j2007T ) with 
the data of lNajita et al.l (|2007l ). The hatched region 
on the plot shows the range of allowed parame- 
ters for a simple p lanet -disc model (adapted from 
Lubow et all Il999l and iLubow fc d'Angelol l2006l ). 
while the grey region marks the range expected for 
"photoevaporated" holes. In the planet-disc model, 
the lines denote the loci where the parameters of 
the disc model are constant: the solid lines are for 
a — 0.01, while the dashed lines denote a — 0.001; 
in both cases the upper line shows the case of a 
IMjup planet, while the low er line denotes a 12Mjiin 
planet. (See discussion in lAlexander fc Armitage 




Fig. 9. Disc masses and accretion rates for discs with in- 
ner holes (adapted from Alexander & Armitage, 200?|). The 
shaded region shows the allowed range of parameters for 
photoevaporative clearing, and the hatched region the range 
expected for a simple model of the planet-disc interaction. 
Th e red and blue poin ts show the transition discs identified 
by iNaiita et al.l ||2007| ): red points are single stars, blue are 
binaries; circles represent detections, triangles upper limits. 
The filled symbols show the objects which satisfy the se- 
lection criterion in Section 4.1; t he open circles show those 
objects classed as transitional bv Na iita et al.l ||2007| ) which 
do not satisfy the more stringent selection criterion out- 
lined here. A l so sho wn are TW Hya (green circle, data from 
Calvet et all l2005h and CS Cha (black cross, data from 
Espaillat et al.Ll2007^ . (Note also that the lower- left of the 
region allowed by the planet model seems rather unlikely 
on formation grounds, as in this region we require that the 
planet be more massive than the disc in which it is embed- 
ded.) 



20071 for more details.) We see that several of the 
canonical transitional discs appear consistent with 
the embedded planet scenario, while only one or two 
of the observed sources fall in t he r ange predicted by 
the photoevaporation scenari However, Fig. 9 



Note, however, that there is considerable disagreement 
as to the significance of the upper limits on accretion 
rates i n cases where no accretion is detected. For ex- 
ample, ISicilia-Aguilar et al.l l l2006bl ) assign upper limits of 
lO~^'^M0yr~^ to objects with no detectable UV excess and 
no broad Ha lines, while the more conservative approach of 
INaiita et al.l (2007) assigns upper limits to such objects in 
the range 10"^''— lO~^M0yr~^. Robust upper limits to the 
accretion rates in objects categorised as WTTs will provide 
valuable additional constraints on models of such systems. 
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also highlights the importance of selection biases 
in analyses of this type. The theoretical arguments 
above consider the formation of deep holes or gaps 
in the inner disc, so it is essential that only objects 
which unambiguously possess such holes be consid- 
ered in the analysis. If we consider the 12 t ransition 
objects identified by iNaiita et al.l (j2007l ). we see 
that 8-10 of them seem consistent with the planet- 
scenario (depending on how one treats the binaries 
in the sample), while only 2/12 are consistent with 
the photoevaporation scenario. If we apply the selec- 
tion criterion suggested in Secti on 4.1, however, we 
reject 4 of the single stars in the lNaiita et al.l ( 2007 ) 
sample, and instead find that 4/6 objects appear 
to be planet-induced holes, while 2/6 are consistent 
with photoevaporation. Moreover, we saw in Sec- 
tion 4.1 that the wavelength range used to identify 
current samples of transitional discs may bias the 
samples towards smaller hole sizes {< 50AU). Given 
that the timescale required for a planet to open a 
gap in a disc is typically ~ 1000 local orbital per i- 
ods (e-g.- IVarniere et al.i . l20o5lEdgar et ani2007bh . 
it seems unlikely that planets at radii beyond a few 
tens of AU will open significant g aps in their discs 
over the Myr lifetimes of TT disc a I . By contrast, 
disc clearing by photoevaporation can produce in- 
ner holes of any size greater than ^ lAU, and while 
the expected hole size distribution depends on the 
viscosity law adopted it is likely biased weakl 
towards larger hole sizes ( Alexander &: Armitagi 
20071 ) . Given the small numbers of transition discs 
known, great care must be taken in applying selec- 
tion criteria to the observed samples, and robust 
samples with uniform selection criteria are essential 
if we are to use such statistical methods to constrain 
mechanisms of disc clearing. 

At present the data on transitional discs do 
not allow us to distinguish clearly between the 
different theoretical scenarios, but these data sug- 
gest that more than one mechanism may pro- 
duce the inner holes we observe. As discussed 



above, the survey of iNaiita et al.1 (|2007l ) suggests a 
planet:photoevaporation hole ratio between 6:1 and 
2:1 (based on very small number statistics). The 
only comparable sam ple of accretion rates i n trans i- 
tional discs is that of ISicilia-Aguilar et al. I (|2006d) . 
who identified 17 transitional discs in the Tr37 clus- 
ter. Unfortunately Tr37 is too distant for TT disc 



masses to be measured, but Sicilia-Aguilar et al] 
(l2006bl ) ibund that roughly half their sample (9/17) 
are accreting, while half (8/17) showed no evidence 
for any on-going accretion. As we observe more 
transitional discs these sample sizes will grow, and 
it seems such studies will provide vastly improved 
statistics within the next few years. 



4.3. Summary of transitional discs 

In this Section we have discussed transitional 
discs, which are thought to be discs observed be- 
tween the CTT and WTT states. We have de- 
scribed the observed properties of these objects by 
considering the best-known sources, and discussed 
the various theoretical mechanisms that have been 
proposed to explain their origin. Infrared observa- 
tions alone do not discriminate strongly between 
different models, and it seems likely that statistical 
studies of disc properties beyond the SEDs offer 
the most promising way of discriminating between 
these models. Lastly, we have discussed various ob- 
servational criteria used to identify these objects, 
including proposing a new, theoretically-motivated 
criterion, and highlighted the strong influence that 
selection effects can have on statistical surveys of 
transitional discs. 



■^"Additionally, models of planet formation and migration 
sugge st that most gas-g iant planets form at radii ~ 5— lOAU 



(e.g., [Armitagi, l2007bl) . 
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